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Abstract 



At the International Linear Collider (ILC), parameters of supersyninictry (SUSY) can 
be determined with an experimental accuracy matching the precision of ncxt-to-leading 
order (NLO) and higher-order theoretical predictions. Therefore, these contributions 
need to be included in the analysis of the parameters. 

We present a Monte-Carlo event generator for simulating chargino pair production at the 
ILC at next-to-lcading order in the electrowcak couplings. We consider two approaches 
of including photon radiation. A strict fixed-order approach allows for comparison and 
consistency checks with published semianalytic results in the literature. A version with 
soft- and hard-collinear rcsummation of photon radiation, which combines photon re- 
summation with the inclusion of the NLO matrix element for the production process, 
avoids negative event weights, so the program can simulate physical (unweighted) event 
samples. Photons are explicitly generated throughout the range where they can be exper- 
imentally resolved. In addition, it includes further higher-order corrections unaccounted 
for by the fixed-order method. Inspecting the dependence on the cutoffs separating the 
soft and collinear regions, we evaluate the systematic errors due to soft and collinear 
approximations for NLO and higher-order contributions. In the resummation approach, 
the residual uncertainty can be brought down to the per-mil level, coinciding with the 
expected statistical uncertainty at the ILC. We closely investigate the two-photon phase 
space for the resummation method. We present results for cross sections and event 
generation for both approaches. 



Zusammenfassung 



Am Internationalen Linearbeschleuniger (International Linear Collider, ILC) konnen die 

Parameter supersymmetrischcr Thcoricn (SUSY) mit cincr cxpcrimentellen Genauigkeit 
gemessen werden, die der Prazision von theoretischen Vorhersagen von nachstfiihrenden 
(next-to-leading, NLO) und hoheren Ordnungen entspricht. Daher miissen diese Beitrage 
in die Analyse der Paramter eingeschlossen werden. 

Wir prasentiercn cinen Monte Carlo Ercignis Generator fiir die Simulation von Chargino- 
Paarproduktion am ILC in NLO in den elektroschwachen Kopplungen. Wir betrachten 
zwei Ansatze fiir den Einschluss von Photon-Abstrahlung. Ein strikter Ansatz fes- 
ter Ordnung (fixed-order) erlaubt den Vergleich und Konsistenztests mit publizierten 
semianalytischen Ergebnissen in der Literatur. Eine Version mit weicher und barter 
kollinearcr Resummation von Photonabstrahlung , weicher die Resummation von Pho- 
tonen mit dem Einschluss des NLO Matrix Elements fiir den Produktionsprozess kom- 
biniert, vermeidet Ereignisse mit negativem Gewicht, sodass das Programm physikalische 
(ungewichtetete) Ereignissamples simulieren kann. Photonen werden in dem Bereich, in 
dem sie experimentell aiifgelost werden konnen, explizicrt generiert. Ausserdem enthalt 
die Methode weitere Korrekturen hoherer Ordnung, die in der fixed-order Methode nicht 
eingeschlossen sind. Durch die Uberpriifung der Abhangigkeit von den Cutoffs, die den 
weichen und den kollinearen Bereich abtrennen, evaluieren wir die systematischen Fehler, 
die infolge der weichen und kollinearen Nahcrung auftreten, fiir NLO Betrage sowie 
Beitrage hoherer Ordnungen. Im Resummationsansatz kann die restliche Ungenauigkeit 
auf Promille-Niveau reduziert werden, welches der experimentellen statistischen Unge- 
nauigkeit am ILC entspricht. Wir untersuchen den zwei-Photon Phasenraum der Re- 
summationsmethodc genau. Wir prasentieren Ergebnisse fiir Wirkungsquerschnitte und 
Ereignisgeneration fiir beide Ansatze. 
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1 Introduction 



In this section, we will give a short overview of supersymmetry (SUSY) and its minimal 
realization, expectations, and results from experiments for SUSY particles, and the 
available computer tools. We also present an outline of the structure of this thesis. 

1.1 Standard Model (SM) and minimal supersymmetric 
extension 

The Standard Model 

The Standard Model (SM) of particle physics [H El El H E] is a SU{3)c x SU{2) x ^7(1) 
gauge theory with the particle content listed in Table [LTl The SU{2) x U{1) part of the 
theory is spontaneously broken by the nonzero vacuum expectation value of the Higgs 
boson [6l [3 O [9l [TO] which gives masses to three of the SU{2) x U{1) gauge bosons. It 
describes current experimental data with high accuracy [11) . However, it suffers from a 
number of theoretical drawbacks. In general, the Standard Model can only be an effec- 
tive low-energy theory as it does not describe gravity and is therefore valid at most up 
to the Planck scale. In addition, it suffers from the fine-tuning or naturalness problem: 
in the Standard Model, corrections to the mass of the Higgs boson are quadratically 
divergent. They can be regularized by a finite cutoff-parameter which can maximally 
be set to the Planck scale as the highest scale of the theory. However, the Higgs mass 
is theoretically bounded [111 |13l HH [151 [16] . Therefore, extreme finetuning is needed to 
obtain the physical Higgs mass from the bare mass in the theory [T7j. Similarly, the 
Standard Model alone does not provide an explanation for dark matter. Prom a more 
aesthetic point of view, it also does not allow for gauge coupling unification of all gauge 
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groups [18]. Further problems are the origin of particle masses, possible symmetries 
between the leptonic and the quark sector, and so on. 

Supersymmetry 

Supersymmetric theories |19l [20] are one of the most promising candidates for the de- 
scription of physics beyond the SM. They introduce a new symmetry as an extension 
of the Poincare group which connects fermionic and bosonic degrees of freedom. This 
symmetry transforms the bosons/ fermions of the Standard Model into their fermionic/ 
bosonic superpartner with a different spin but otherwise identical quantum numbers. 
As the superpartners have not been directly observed in experiments, SUSY has to be 
broken such that all new particles obtain higher masses. The breaking most likely takes 
place in a "hidden sector" invisible to the Standard Model gauge groups and is then 
transferred to the visible sector. There are numerous suggestions for SUSY breaking 
mechanisms, including (minimal) supergravity [211 [22l [23] , gauge-mediation [2 H [25 [ |26]. 
and anomaly-mediation |27[ [28] . An additional symmetry (R-parity) only allows for the 
simultaneous creation/ annihilation of even numbers of SUSY particles. R-parity viola- 
tion is strongly limited by experimental results from proton decay. 
Supersymmetric theories address two main problems of the Standard Model: they allow 
for a natural solution for the fine-tuning problem. In addition, if R-parity is exactly 
conserved, the lightest supersymmetric particle (LSP) is a good dark matter candidate. 
Furthermore, they easily allow for the embedding of the SU{3)c x SU{2) x U{1) group 
in a grand unified theory and the unification of gauge couplings. 

For more technical details on the construction of supersymmetric theories, cf. Appendix 

El 

Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric Standard Model (MSSM) is the minimal supersymmetric 
extension of the Standard Model |29| , 1221 [30l 131] . It contains a superpartner for each SM 
particle and an extended Higgs sector with the two Higgs doublets = (^o ) , = 

[it-)- The Lagrangian at the electroweak scale is given by 

d 

^tot — ^kin ~l" ^ gauge ~\~ ~l~ C^softt 

where Ckin contains the kinetic terms of the free theory, Cgauge are terms arising when 
imposing the gauge symmetry of the SM, Cy is the Lagrangian part derived from the 
superpotential, and Cgoft are the soft SUSY breaking terms. For more details and the 
complete MSSM Lagrangian, cf. Appendix \K\ While Ckm and Cgauge only depend on 
SM parameters, new SUSY related parameters appear in the superpotential and the 
soft breaking terms. The unconstrained model contains in total 105 new parameters. 
The number of free parameters, however, can be constrained by imposing lepton number 
conservation, suppression of flavor changing neutral currents, and applying experimental 



5 



bounds on CP violation. The assumption of a specific breaking mechanism can reduce 
the number of new parameters even further. In the mSugra scenario of the MSSM, all 
parameters can be derived from 5 new parameters at the SUSY breaking scale. The 
parameters at the electroweak scale are then determined using renormalization group 
equations. 

For recent reviews, see [5 ^ [55 1 [M]. 

The Chargino and Neutralino Sector of the MSSM 

A solid prediction of the MSSM is the existence of charginos xf:X2 neutralinos 
X?) X2' X4- These are the superpartners of the and charged/ neutral Higgs 
bosons H^, rotated into their mass eigenstates. In grand unified theories (GUTs) the 
chargino masses tend to be near the lower edge of the superpartner spectrum, since the 
absence of strong interactions precludes large positive renormalizations of their effective 
masses. The precise measurement of chargino parameters (masses, mixing of with X2, 
and couplings) is a key for uncovering any of the fundamental properties of the MSSM. 
These values give a handle for verifying supersymmetry in the Higgs and gauge-boson 
sector and thus the cancellation of power divergences. Charginos decay either directly or 
via short cascades into the LSP, which is the dark-matter candidate of the MSSM in case 
of R-parity conservation. Thus, a precise knowledge of masses and mixing parameters 
in the chargino/neutralino sector is the most important ingredient for predicting the 
dark-matter content of the universe. Finally, all SUSY parameters in the gaugino and 
higgsino sector of the MSSM can be determined from the measurements of the chargino 
production cross sections and masses and the lightest neutralino mass |35j . The high- 
scale evolution of the mass parameters should point to a particular supersymmetry- 
breaking scenario, if the context of a GUT model is assumed (cf. [MIEZ]). In all these 
cases, a knowledge of parameters with at least percent-level accuracy is necessary. 

1.2 SUSY at colliders: discovery and precision 

In this section, we give a quick overview on SUSY searches at past and future colliders. 
We refer to [38] and [391 SO] and references therein for details. 

Experimental bounds 

Direct experimental searches for SUSY particles have been conducted at both LEP and 
the Tevatron. Neither experiment can claim a direct discovery for a SUSY particle, but 
both have determined lower mass limits ^38j. While the combined LEP runs provide the 
most stringent mass limits for sleptons and gauginos, limits for squarks and gluinos can 
be obtained from CDF at the Tevatron and ALEPH at LEP. Most lower limits for visible 
SUSY particles (ie, not the LSP) are 0(100 GeV). However, limits also depend strongly 
on the underlying SUSY scenarios and breaking mechanism. For a recent review for the 
mSugra parameter space, cf. [H]. Results from HERA give limits on R-parity violating 
scenarios [32 l Ii3l liS] . 
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Indirect constraints on the SUSY parameter space are in general given by (loop-induced) 
effects in electroweak precision data. The most stringent restrictions arise from the mass 
of the lightest Higgs boson. Further observables include the W boson mass, sin^ 6w, the 
flavour changing neutral current 6 — > 57, and the anomalous magnetic moment of the 
muon {g^ — 2) [35]. 

Additional restrictions on the SUSY parameter space can be obtained from cosmological 
data, as dark matter searches. 

LHC: discovery 

The LHC as a hadron-hadron collider with a cm energy of ^/s = 14TeV serves as a 
discovery machine for SUSY particles in most SUSY scenarios. The most dominant 
channels are pair production or associated production of gluinos and squarks. As su- 
persymmetric particles usually decay via long decay chains, mass determination involves 
reconstruction from combined mass distributions. Analyses are parameter-point specific. 
For the mSugra point SPSla(')[l6](cf. Appendix [C]), analyses for chargino, neutralino, 
squark, gluino, and slepton masses have been done \39\ I47j. Relative errors for masses 
and cross sections are usually expected to lie in the % regime. 

International Linear Collider (ILC): precision 

The International Linear Collider (ILC) is a planned e"*" e~ machine with a cm energy 
of 500 GeV (ITeV). It provides much cleaner production channels and decay signatures 
with lower background than the LHC. In addition, direct pair production of sparticles 
in the slepton and chargino-/ neutralino sector is easily observable. If kinematically 
accessible, the sparticle content of the MSSM can be studied with high precision [38| 
I39j . Combining sparticle measurements from the ILC and the LHC leads to errors 
in the percent to per mille regime and significantly reduces the errors based on LHC 
measurements only [491 I39| . The same accuracy is reached in the respective fitting 
routines Sfitter [50] and Fittino [51] for the reconstruction of the Lagrangian parameters 
at the weak scale. 

Due to its low mass in most SUSY scenarios, the lighter chargino is likely to be pair- 
produced with a sizable cross section at a first-phase ILC with cm. energy of 500 GeV. 
In many models, including supergravity-inspired scenarios such as SPSla/SPSla' [16], 
the second chargino will also be accessible at the ILC, at least if the cm. energy is 
increased to 1 TeV. Similar arguments hold for the neutralinos. 

1.3 Next-to-leading order at Monte Carlo Generators 

Several event generators [52l [531 [Ml [55] already include the particle content of the MSSM 
and allow for event generation including sparticles at tree level. In Ref. [56], some of 
these have been presented and verified against each other, both for the SM and the 
MSSM. 
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To match the experimental accuracy for (SUSY) processes at (future) colhders, theoret- 
ical predictions with next-to-leading order (NLO) accuracy have to be implemented in 
the simulation tools (see e.g. |40j ) . The inclusion of higher-order corrections in Monte 
Carlo event generators has already been discussed for electroweak precision physics at 
LEP (cf. \57 \ [58 l [59]). For the determination of the W-mass, Monte-Carlo event gener- 
ators including higher-order corrections \tjO\ [UT] can reduce the theoretical uncertainties 
down to approximately 0.5% [62 | [63 | ITT]. 

NLO calculations for the production of charginos and neutralinos at an e'^e' collider are 
within the percent regime [64 ^ 1651 [66] . If SUSY is discovered, the experimental precision 
at the ILC requires to equally include these processes at next-to-leading order in a Monte 
Carlo event generator. Furthermore, it is essential for the simulation of physical (i.e., 
unweighted) event samples that the effective matrix elements are positive semidefinite 
over the whole accessible phase space. The QED part of radiative corrections does not 
meet this requirement in some phase space regions. This usually requires higher-order 
resummation for soft (and collinear) photonic contributions to the process. Methods for 
dealing with this problem have been developed in the LEPl era |67 [ I68 [ [57]. Since the 
ILC precision actually exceeds the one achieved in LEP experiments, these higher-order 
effects from resummation can have the same order of magnitude as the experimental 
errors. 

Similar techniques concern the event simulation of hadronic processes including partonic 
showers. The higher-order parton showers have to be matched with the exact NLO ma- 
trix element. Work along these lines has been done for e~^e~ |69ll70j and hadron colliders 

[ZI11Z211Z3]. 
1.4 Outline 

In this work, which is the extension of a recent publication [71] , we present an extension 
of the Monte Carlo Event generator WHIZARD [52J which includes the 0{a) electroweak 
and SUSY corrections to chargino production at an e"*" e~ collider. The fixed-order ver- 
sion, which we will discuss first, is limited to first order photon emission and suffers the 
well-known problem of negative weights for low cuts on the photon energy. We then 
introduce a method which combines the photon resummation already discussed at LEP 
with NLO matrix elements for the production process. It avoids negative events and 
therefore allows for lower energy cuts. In addition, it includes further higher-order cor- 
rections unaccounted for by the fixed-order method. We evaluate the systematic errors 
due to soft and collinear photon approximations and the cut-dependence of higher-order 
contributions. We present results for the cross section and event simulation for both 
methods. 

In Chapter [2], we discuss chargino production at an e"*" e~ collider at tree level, pre- 
senting the total and differential cross sections for any helicity combination. In Chapter 
m we present the analytic form of the infrared and cut-independent cross section for 
chargino production at next-to-leading order and the inclusion of resummed soft and 
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virtual photonic contributions. We sketch the treatment of 2 — > re processes including 
chargino decay in the double pole approximation. 

In Chapter m we present a method to include the fixed 0{a) result in the Monte Carlo 
event generator WHIZARD. We give a brief overview on the technicalities of the inclusion 
and discuss the energy-cut-dependent problem of negative weights in certain points of 
phase space. We present first results for cutoff-independence, cross section integration 
and event generation in the allowed cut parameter regions. 

In Chapter O we present our method of combining the completely resummed soft and 
virtual photonic contributions with the exact next-to-leading order contributions to 
chargino production. We explicitly discuss the description of collinear photonic con- 
tributions up to 0(0^) and analytic differences between the resummation and the fixed 
order method. We present the results of our work in Chapter [6l In Chapter [TJ we sum- 
marize and give an outlook on future work. 

The appendices basically contain conventions, derivations of SUSY Feynman rules and 
the approximations used for the photonic contributions. We start with a general dis- 
cussion of supersymmetric theories and the MSSM in Appendix El The chargino- and 
neutralino-sector, especially the derivation of Feynman rules, is discussed in Appendix 
IBI Appendix|D] gives an overview on the treatment of helicity states for massive particles 
needed for the discussion of the helicity-dependent Born cross section. 
Appendix [El shows all generic diagrams contributing to chargino production at NLO. 
Appendix [Fl sketches the derivation of the collinear photon approximation and the ini- 
tial state radiation (ISR) structure function [75]. We present the mSugra point SPSla' 
in Appendix O and list some useful functions in Appendix [Gl Appendix |Hl lists the 
references for all computer programs mentioned in this work. 
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2 Chargino production at tree level 



In the MSSM, charginos x ^-re superpositions of superpartners of the charged gauge 
bosons and the charged Higgs fields H^. In the fohowing, we use the mass- 
diagonahzation matrices 



\ -sin$L,R cos ^L,R J 
for the case with no CP violation. For more details, cf. Appendix iBl 

2.1 Born matrix element 

At an e~*" e~ collider, charginos are produced by 7 and Z exchange in the s-channel and 
z>e exchange in the t-channel as shown in Figure 12.11 Neglecting the electron mass nie, 
we can perform a Fierz transformation on the t-channel De exchange. The sum of all 
contributions is then given by [76] : 



with a, (3 = L,R and P^//? = ^ {I ^ 75). The factors Qafj have been explicitly 
calculated in ^76j and are given by 




and the parameterization 




(2.1) 



Qll = DlT Fl cos 2$l 
Qlr = D'lT F'l cos 2^r 



Qrl 
Qrr 



Dr =F Fb^ cos 2^l, 
Dr =F Fr cos 2 ^R, 



X 



,+ 



X 



e 



e 




e 



X 



X 



Figure 2.1: Feynman graphs for chargino pair production at the ILC. 
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for the production of xf Xi ^-nd xt X2 

Qll = Fl sin 2^>l , Qrl = Fr sin 2$^, 
Qlr = F'l sin 2 , Qrr = Fr sin 2 $r 

for the production of xf X2 xt Xi pairs. We here use the notations [75] 



Dr 



F'l 



1 + 



1 



,2 



Dz 



Fl-^, Dz 



4s 



w 



rriz + tTz 



, Dr. 



Dz 



4,2 2 



4,2 ■ 



m 



2 ' 



where s^, and denote the squared sine and cosine of the weak mixing angle. 



(2.2) 



2.2 Helicity amplitudes 

The matrix element Ai given by Eq. (|2.1|) can be decomposed into helicity amplitudes. 
We use the formalism introduced in [77] (see Appendix [D] for a short review on helicity 
eigenstates and an introduction to the used method). We define the direction as the 
flight direction of the electron and obtain 



Pi(e+) 



P3(X" 






V -Pe J 

Ex \ 
sin 6 



V Py^cosO j 



P2{e ) 



P4(x" 



/ Pe \ 




\Pe / 

( Ex 
—p^ sin 


V -p^cosO J 



\ 



(2.3) 



where pe and py. are the magnitude of the electron/ positron and chargino three- momentum, 
respectively. 9 is the angle between the electron and the negatively charged chargino. 
M{a; Xi Xj) = 2 n a {a; Xi Xj) defines the matrix element for an electron with helicity a 
and with helicity Aj/j. The helicity amplitudes are then given by [78] : 



(+;++> 

(+;+-) 
(+;-+) 

(+;— > 



Qrr \n-rii + Qrl \P~-V 



smi 



Qrr ^(1 + r]+) (1 + ??-) + Qrl V(1 " V+) (1 - J (l + cos0), 

'^ + QRLVi^ + V+){l+V-)] (1 



QrrV(1 
Qrr 



1 



cos ( 



i-vi + QRL \n-v 



smt 
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(- 


1 1 \ 

++) 


\- 




(- 


-+) 


(- 


— ) 



sin( 



Qlr V(1 + ^+)(1+»7-) + Qll (1 - cos 9) , 

Qlr V(1 - r/+) (1 - + Qll V + ri+) {I + V-)] (1 + cos 9) , 



QLRjl-rii + QLLjl-lli 



sin( 



with 



r]± = VX ± {i4 - nj); \ = [l - {m + Hjf] [l - - Hjf] ; ^ 

2.3 DifTerential and total cross sections 



For the calculation of the differential cross section, we use the general formula for 2 — 2 
particle processes: 



da 



1 d^ki 1 



2w 2k1 (27r) 



S{kl-m'i)Q{k^^)\M\\ 



where ki, k^, are the momenta/ energy/ mass of the outgoing particles. The flux 
factor 

w{s,ml,ml) = -y/^ 



il)^ -4ml 



( s — m| — m 

depends on the center of mass energy -y/s and the masses of the incoming particles. For 
particles with very small masses with respect to the cm energy, w = s. 
The differential cross section for a specific helicity state is then given by 



da 



TT a 



VX\{a;\iXj)\' 



d cos 9 8 s 

For chargino pair production, the differential cross section is 
da --'2 



dcos 6 
with 



^ \/A { [1 - {fif - ii]f + A cos^ 9] Qi + AnifijQ2 + 2 cos 9^ (2.4) 



Qi = I [\QRRf + IQllI' + \QRLf + \QLRf] , 
Q2 = ^Re [Qrr Q\l + Qll QIr] , 

Qs = \ [\QRRf + \QLLf - IQrL^ - \QLRf] , 
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-1 -0.5 0.5 1 

cos 9 



Figure 2.2: Chargino pair production at the ILC: Dependence of the differential distri- 
bution in polar angle cos 9 between e~ and Xi for different helicity combina- 
tions. The labels indicate Xi xt helicity; the electron/positron helicity 
is fixed to — h. 

where we averaged over initial and summed over final particle spins. Figure 12.21 shows 
the results for the differential cross section for the dominant helicity amplitudes for 
xt Xi production for the point SPSla' (cf. Appendix [C|). Independently of the SUSY 
parameters, processes with a left handed electron in the initial state are dominant, as 
diagrams with a 7 and Z exchange interfere constructively. The same diagrams give 
destructive interference between the 7, Z exchange for righthanded initial electrons. 
This can easily be seen from Eqs. (|2.2|) . We obtain 

Dl ~ I + 0.727 Dz, Dr ^ I- 0.67 Dz, Fl « -0.01 Dz, Fr k, O.SSDz- 

For the cm values where chargino pair production is possible, Dz > 1. Then, even 
without a possible enhancement from the t-channel sneutrino exchange and independent 
of the mixing parameters, \Qrl/Qll\ < 0.33. The ratio of \Qrr/Qlr\ strongly de- 
pends on the mixing parameters and the sneutrino mass. For the point SPSla' and a 
cm energy of 1 TeV (400 GeV), \Qrl/Qll\ = 3%(1%) and \Qrl/Qll\ = l%o(l%) (ne- 
glecting contributions from the sneutrino exchange). We therefore only show the results 
for states with a left-handed electron in the initial state. Amplitudes with two charginos 

of the same helicity in the final state are suppressed by a factor oc . Figure 12.31 shows 
the helicity averaged/ summed differential cross section according to Eq. ()2.4p . Figure 
12.41 the ^/s dependence of the corresponding total cross section. 

For the same parameter point, we have compared the tree level results for Xi Xj Pro- 
duction using the analytic form [76j and several computer codes [52l \79\ [HOl [81] • The 
analytic and numeric results are in complete agreement. The values are given in Table 
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cos 6 

Figure 2.3: Chargino pair production at the ILC: Dependence of the spin summed and 
averaged differential distribution in polar angle cos 9 between e~ and Xi ■ 
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rr, , \ fh] 


Xi Xi 


ioD.O 


XiXt 


10.97 


X2Xt 


80.07 



Table 2.1: Chargino total production cross sections for the SUSY parameter set SPSla' 
and = 1 TeV. 

2.4 Polarized incoming and outgoing particles 

In |76tl78j. the authors also list the analytic results for the x~ polarization vector as well 
as the e^-polarization dependent differential cross sections. Although not pursued fur- 
ther in this work, we reproduced these results and sketch the derivation for completeness. 



Polarized outgoing particles 

The polarization vector of the outgoing x~ iii its rest-frame is defined in the following 
coordinate system: z is given by the component parallel to the flight-direction of the 
chargino in the lab-frame, x is defined in the production plane, and y normal to the 
production plane. The x~ polarization vector is then given by 

V =tr{-^p) = {Vr,VM,Vc), 
where ~a is the vector made of the Pauli- matrices Cj , p the spin-density matrix, and 

Vr = ^i?e|5](a;++)(a;-+)* + (a;— )(a;+-)*|/Ar, 
Vm = )(a;+-r - (a;++)(a;-+)*|/AA, 

^ = iE [K+;A.A,)p + |(-;A.A,)n. 

\i Xj 

Polarized incoming particles 

For the calculation of the differential and total cross sections depending on the polariza- 
tion of the incoming particle beams, we start with 

V = {Vr, 0, Vc) , V = {V'r cos r?, V'r sin r?, -V'c) 

for the polarization vectors of the electron and positron, respectively, in a coordinate 
system where the z-axis is given by the momentum of the electron and the x-axis by 
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the electron's transverse polarization vector, rj is the azimuthal angle of the positron 
transverse polarization vector with respect to the x-axis. We transform this into the 
lab-system, where the x-axis is defined by the scattering plane (cf. Eqs (j2.3p ): 

r = {Vr cos ^,-'Pr sin <^,'Pc) ; V = (Pr cos(r/ - ^>),P^sin(?7 - ^),-Vc); 

$ is the angle between the transverse polarization vector of the electron and the new 
X-axis. 

The squared matrix element including an incoming particle with possible helicity eigen- 
states labeled by A is given by [82] 

X,X' 

for particles and 



A,A' 



for antiparticles. The spin density matrix is given hy p = ^{l + 'aV), and V is defined 
such that J) = \p\ez- Taking this into account, the electron spin density matrix is 

1 / l + Vc Vre'"^ 



For the positron, we still have to perform a rotation around the x-axis such that p'^ 
\Pe+\ez,e+ leading to 

P' = (p^cos(T/-$),-P^sin(r/-$),P^) 

which then gives 



If we now calculate the sum over the squared helicity amplitudes using 



a'i,(72,(TjO-2;Ai 

and 

M{aia'i;Xj) = M{ai; Xj)6„^ 

we obtain 
da G? 



dn 16 s 
with 



VX [(1 - rcr'c)^unpoi + {Vc - V'c) + Vr V'r cos (2 $ - r?) Stt] 



-'unpol 



4 I [1 - (/uf - p]f + A cos^ e] Qi + Api pj Q2 + 2VXQ3 cos^} 
= 4^^^[l-{p^-p]f + Xcos^d] Q[+ 4 Pi pjQ'^ + 2 VXQ's cos 9^ 
T,TT = -4Xsm^eQ5, 
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Q'l = \ [\Qrr? - \Qll? + \Qrl? - \Qlr?] , 

Q'2 = 2 [Qrr Q*rl - Qll Qlr] > 

Q3 = \ [\Qrr\^ - \Qll\^ - \Qrl\^ + \Qlr\^] , 

Qb = ^Re [Qlr Q*rr + Qll Qlr] 

and Hi J and A as introduced in the last section. In [76l [78], the authors show that the 
measurement of the hghtest chargino mass, the production cross sections, and spin-spin 
correlations suffice to completely determine the parameters of the chargino system in 
the tree level approximation. 
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3 Chargino production at next-to-leading 
order (NLO) 



3.1 Divergencies and infrared-safe cross sections in 
higher-order calculations 

Before discussing the 0{a) corrections to chargino production at an e~ colhder, we 
Ust some general features of the calculation of (virtual) higher-order corrections in per- 
turbation theory. 

In finite order calculations, divergences in the ultraviolet {E oo) as well as infrared 
(E — > 0) regime can appear. The UV divergencies have to be regularized, which in- 
troduces a regularization parameter A. This leads to a parameter-dependence of the 
physical quantities with respect to the bare parameters in the Lagrangian. In bare per- 
turbation theory, the bare parameters are then eliminated in relations between physically 
measurable quantities. If the underlying theory is meaningful, these relations should be 
cutoff-independent. Alternatively, the theory can be renormalized by absorbing the UV 
divergencies in a redefinition of the parameters and fields of the theory, thus giving a 
physical meaning to the parameters in the Lagrangian. The bare Lagrangian is split into 
the renormalized Lagrangian and a counterterm part, 



where the Feynman rules resulting from Cct have to be included in all calculations of 
physical observables. 

IR divergences arise in the case of zero-mass virtual gauge boson exchange. When virtual 
massive gauge bosons with the mass rUg are connected to an on-shell particle with the 
mass m, logarithms of the form 



The IR divergence in QED can be regularized by introducing an infinitesimal gauge 
boson mass A. According to the Bloch-Nordsieck theorem [83j, the divergencies cancel 
if the emission of soft real photons is also taken into account. Infrared-safe observ- 
ables therefore always include terms describing the emission of soft real photons. In 
the collinear approximation, where the transverse momentum of the emitted photon is 
neglected, the dominant contributions originating from the multiple emission of soft and 
virtual photons from the initial particles can be summed up into structure functions. 
This is discussed in Section [3771 



■bare — C^ren ~l~ l~'cti 




appear. This term gets infinite as m, 



'9 



0. 
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In addition, collinear divergencies can appear when the (massless) gauge boson is emit- 
ted at a small angle relative to the emitting particle. In the collinear approximation the 
integration over the photon phase space leads to terms of the form 

,max 

log 

m 

which diverge when m = 0. These singularities are regularized by keeping the physical 
nonzero mass m in this region of phase space, cf. Section [3.41 

An infrared-safe total cross section with a cm energy y/s and n particles in the final 
state then includes the following contributions 



Born cross section: 



(TBorn{s) = j dTn\MBorn?-, (3.1) 



• interference term between Born and first order terms describing the purely virtual 
contribution 

CTvirt(s, A2) = j dVn [2Re {MboUs)* A^l-loop(s, A^))] (3.2) 

• soft photon contribution 

Born (^)r (3.3) 

(The soft photonic approximation is explained in Section [3.3p . 

Here, A denotes the photon mass and A is the soft photon energy cutoff separating 
the hard from the soft region, fsoft (!3.6p denotes the soft photon factor. 



The sum of the three contributions ()3.ip . ()3.2h . ()3.3p 

a'{s,/S.E^) = aBorn{s) + (Tvirt(s,A^) -hfTs(s, AE'^, A) 

is infrared-safe. However, it still depends on the soft photon cutoff I\E^ as the soft 
approximation only takes a part of the m ^ n + 7 phase space into account. For a 
cutoff-independent result, the hard cross section 



— * (n + 1) 

Jae^ 



2 



has to be added. This is usually split into a hard, collinear and a hard, non-collinear 
part 

<7m ^ n + l(s, AE^) = (Thard,non-coll(s, A_E^, A6'^) + (Thard,coll(s, A£'-y, A^^), (3.4) 
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where the cutoff A 0^ separates the colhnear from the non-colhnear region. The hard- 
colhnear part is treated in Section f3.4l the non-cohinear in Section [3.51 
The total 0{a) cross section 

crtot{s) = cFBorn{s) + cr^{s,\) + fJs(s, A£'^, A) + (Jm^(„+i)(s, A£;^, A6'^) (3.5) 

is then cutoff- independent. 

3.2 Virtual corrections 

The one-loop corrections to the process e~ XiXj with i,j = 1,2 have been com- 
puted in the SUSY on-shell scheme in Ref. [641 . An independent calculation in the 
DR scheme has been presented in [Sn]. These calculations include the complete set of 
virtual diagrams contributing to the process with both SM and SUSY particles in the 
loop. The collinear singularity for photon radiation off the incoming electron/positron 
is regulated by the finite electron mass mg. As an infrared regulator, the calculation 
introduces a fictitious photon mass A. Both calculations use the FeynArts/FormCalc 
package [ 85^ [HOj ISTl [86] for the evaluation of one-loop Feynman diagrams in the MSSM. 
A complete list of all generic one-loop diagrams (excluding tadpoles) contributing to 
the process e"*" e~ xt Xi is given in Appendix |El The total number of diagrams is 
0(1500), with 0(1000) self-energy diagrams for all contributing particles. It also in- 
cludes e'''e~Higgs couplings which are absent if rn-g is set to zero. We refer to [641 184^ [66] 
for details of the calculation. 

3.3 Soft terms 

The soft-photon factor has been derived in [HTl [88l [89l [90] . We just sketch the derivation 
and refer to the literature for further details. 

In the soft photon approximation, the radiation of a photon off an incoming or outgoing 
charged particle for an arbitrary process is considered in the limit — > 0, where only 
the terms contributing to the infrared singularity are kept. Then, the radiation can be 
described by a factor depending on the particle and photon momenta and the matrix 
element describing the radiation is proportional to the Born matrix element: 

where Pe± (^) denote the electron/positron (photon) four- vectors. For an incoming 
fermion, this factor is given by 

E p 

f{p,k) = -^Qf-j^- 

Here, Qi? is the charge of the fermion and e the polarization vector of the photon. In 
addition, 

ujk = Vk^ + A2 
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is the energy of a photon regularized by the photon mass A. 
For cross sections, this leads to 



(^m — > n- 



,+7 — f softer-, 



m 



where 



/soft 



27r 



a 



E 



f (fk^ {±)ptPjQiQj 

|k|<AE^ 2u;fc {Pik){pjk) 



(3.6) 



is summed over all charged incoming/ outgoing particles. The (it) in the numerator 
depends on the charge flow (— for an incoming and + for an outgoing charge). The 
integral appearing in fsoft has been calculated in [RTj and is given in Appendix IF.ll 
asoft is then given by 



3.4 Hard-collinear photons 

If a photon is emitted off a particle of mass m with a small transverse momentum k±, 
logarithms of the form 



appear. Considering photon emission off electrons, these logarithms give rise to diver- 
gences if the electron mass is set to zero. Therefore, in order to regulate these collinear 
divergencies, the electron mass has to be taken into account in these regions of phase 
space. Similarly, numerical integration becomes tedious or even unreliable for very small 
transverse photon momenta. It is therefore customary to use an analytic collinear ap- 
proximation for small k± (or, alternatively, small emission angles 6^) when integrating 
over these regions of the photon phase space. 

The explicit expression for the (hard) collinear photon approximation has been derived 
in [9I1E21I93]; cf. also Appendix El 

The hard-collinear contribution to the cross sections from the radiation of photons off 
one incoming particle is then given by convoluting the Born cross section with the struc- 

2 

ture function f'^{x; A9^,^), with x = 1 — 2E^I^/s being the energy fraction of the 
electron after radiation. 




(3.7) 






(3.9) 
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e+ X e 



Figure 3.1: Matrix elements contributing to the 2^3 process. We only show the s 
channel contribution (t channel exchange analogously). 



The two structure functions /"*" , / are given by (jF.Sp 

^ ^ > 27r (l-x) V V J J ' 

Z vr 

They correspond to helicity conservation and helicity flip, respectively; each one is 
convoluted with the corresponding matrix element. The cutoff AE^ is replaced by 
xo = 1 — 2 AE^/^/s. In this approximation, positive powers of A9^ are neglected. For 
radiation off both incoming particles, we have 

0"hard, coll = <7hard, coll(e''') + CThard, coll(e )• 

3.5 Hard non-collinear photons 

The hard non-collinear contributions are added in form of the analytic e"*" e~ x~ 7 

matrix element. In order to prevent double counting for soft and hard-collinear photons, 
which are already accounted for in Osojt and Ohard.coii , lower angular and energy cuts for 
the explicitly generated photon are set. Then, 

(^hard,non-coll = / dPs |7W2^3 (•?) P • (3.10) 
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In the non-collinear region, logarithms as (j3.8p . which are regulated by a finite photon 
mass, do no longer appear. For cm energies 0(100 GeV) and larger, we can neglect 
contributions proportional to the electron mass and set mg = 0. 

The contributing Feynman diagrams for the s channel exchange are shown in Figure 13.11 
The matrix elements can be easily obtained from Eq. (12. ip by substituting 

M = Au(p) — > M' = ^ A{i)-}^ + m)iu{p), 

2pk ^ 

M = A'v{p') M' = --^A'{]^' -]^ + m)^vip'), 

Ip'k 

where u or v are the spinors of the particle radiating off the photon with a momentum 
and the polarization vector and A /A' symbolize the part of the matrix element 
untouched by the radiation. 
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Figure 3.3: relative difference between Born and fixed-order result for total cross section: 

^ Born 



3.6 Fixed 0{a) results for total cross section 

The total fixed-order 0{a) cross section (Eq. (13. 5p ) for the process e~ —>■ xt Xi ^ (t) 
is then given by the sum of Eqs. ([31]), (1^21), <^ 

Crtotis) = O-Born(s) + CTsoftis) + Crvirt(s, , m^) + (Thard,coll(s, AE^^, A6'^, m^) 

+ O'hard,non-coll('5> AE'^, A6'^) 

= j dr2\MBoTn{s,COSe)\'^ + j /soft |XBorn(s, COS 0) 

+ j [2Re {MBornisT A^l-loop(s, A^, ml))] 

+ V V r dx,r{xi;Ae^,^) [ dT2\M^,,,{xi,s,ml)\^ 

+ [ dT3\M2^3{s)\'' 

(3.11) 

with the contributions introduced in the previous sections. For the mSugra point SPSla' 
(cf. Appendix [C]), this leads to corrections 0(20%) near the threshold and in the 5% 
region for ^/s > 600 GeV; cf. Figures [3^2] and [331 



3.7 Resummation of higher logarithms: Initial state radiation 

The logarithms (|3.8|) which arise in the collinear emission of photons can become large 
for small m. In the collinear approximation, where the transverse momentum of a photon 
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with respect to the emitting particle is neglected, the divergencies originating from the 
emission of real collinear photons as well as their virtual counterparts can be summed up 
in splitting functions. The electron-electron splitting function in the leading logarithmic 
approximation is given by 

Peeix) = (^[^)^' (3-12) 

where only terms proportional to the collinear logarithm are kept. The +-distribution 
is given by Eq. ()F.6p . The differential cross section taking the emission of one real and 
one virtual photon collinear photon into account then reads 



(x) dax{xp), 



where p is the electron momentum, X symbolizes the final state of the reaction and Q 
is the scale of the process. The distribution function 

/e,e(x,Q2) = S{l-x) + ^ log (^) Peeix) 



gives the probability of finding an electron with the longitudinal momentum fraction x 
in an incoming electron, when the emittance of photons with the maximum transverse 
momentum are taken into account. In the collinear limit, Q should be set equal to 
k± and be small compared to po- 

If the coherent emission of more than one photon is considered, the splitting function is 
replaced by a distribution function D^^ . The dominant logarithmic contributions stem 
from the emission of photons with strong A;_L-ordering such that k±^i ^ k±^2 ^ •••• 
The distribution function D'^^ then obeys the evolution equation 



|.0-,.„,.^£^P(.,0«.(£„,, (3,13) 



where 

rj 



ds'2a{s') 



s vr 

In a first approximation, we can neglect the running of a and set it constant. In this 
case, 

^ = — log (^) . (3.14) 
The first order solution of Eq. (j3.13|) is then 

I)^^(x,r/) = 5{l-x) + '^Pe,e{x) + 0{lf), 

where the (5(1 — x) term corresponds to the tree level (= no photon emission) part, 
higher-order solutions can be found in the literature; cf. \Jb \ 195^ [96]. The exponentiated 

2 

structure function /isR(a;; A0-y, ^) (Eq. (IF.14P ) includes photon radiation to all orders 
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Figure 3.4: Relative effect of ISR higher-order initial state radiation: 



'^tot.ISR(b)-<^tot 



in the soft regime at leading-logarithmic approximation and, simultaneously, correctly 
describes collinear radiation of up to three photons in the hard regime. It does not 
account for the helicity-flip part /~ of the fixed-order structure function. More details 
on this can be found in Appendix IF.3[ 

Initial state radiation from both incoming particles is then given by 

CrBoTn+lSR{s,^0^,ml) = j dxifiSRixi) j dx2 hsK{x2) j |A4Born(a:i,X2,s)|^. 

(3.15) 

In combining utot (Eq- I3.11|l and UBom+iSR (Eq. I3.15p . we have to subtract the 0(1) 
and 0{a) contributions of OBom+iSR as they are already accounted for by a tot- This 
way, the total cross section including all 0{a) contributions as well as higher-order initial 
state radiation, is given by 

(^tot,ISR{b) = f^tot + (^Born+ISR — Born 

- i^j (ixi /i'^r(xi) j dEa |7WBorn(xi,s)|^ + j dx2fisR{x2) J dT2\MBorn{x2,' 

(3.16) 

where ffgji is the 0{a) contribution of /isr- NLO results for total cross sections in 
the literature are typically presented in the form of Eq. (j3.16p . The effects of including 
higher-order ISR are in the per mille regime as can be seen from Figure! 



3.8 Further higher-order contributions 

Charginos usually decay via decay chains involving (at least) two final state particles. A 
complete NLO calculation therefore also includes factorizable corrections to the chargino 
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Figure 3.5: General example for a non-factorizing NLO contribution to the 2 — > 4 
process with s channel chargino production and unspecified final states 

decays as well as non-factorizable corrections as e.g. in Figure [531 

A first step is the use of the double-pole approximation for the unstable particles. It 
includes (i) loop corrections to the SUSY production and decay processes, (ii) nonfactor- 
izable, but maximally resonant photon exchange between production and decay, (iii) real 
radiation of photons, and (iv) off-shell kinematics for the signal process. Recent com- 
plete NLO calculations for SM W pair production at an e~^e~ collider [97] have explicitly 
verified the validity of this approximation in the signal region. A complete analysis also 
requires the consideration of (v) irreducible background from all other multi-particle 
SUSY processes, and (vi) reducible, but experimentally indistinguishable background 
from Standard Model processes. So far, no calculation provides all NLO pieces for a 
process involving SUSY particles. 

In this work, we only consider the extension of the tree-level simulation of chargino pro- 
duction at the ILC by radiative corrections to the on-shell process, i.e., we consider (i) 
in the above list and consistently include real photon radiation (iii). This is actually a 
useful approximation since nonfactorizable NLO contributions are suppressed by 0{^) 
[98\ 199] and in many MSSM scenarios the widths of charginos, in particular x^, are quite 
narrow (cf. Table [3^1]) . 

Chargino decays (ii), non-factorizing dominant contributions (iii), and finite- width ef- 
fects (iv) will be covered in future work. NLO corrections to chargino decays for specific 
decay products are available from [84] . They are in the % regime and can easily be com- 
bined with the present analysis. The inclusion of background effects (v) and (vi) can 
easily be done using Monte Carlo Event generators [56]. Non-factorizing contributions 
and finite width effects can be treated in the double-pole approximation |10Ul IIUI] [BO]. 
Here the propagators of the unstable particles are expanded around their poles, and only 
leading order contributions are kept. For W pair production at an e^e~ collider, ana- 
lytic results for non-factorizing contributions [1021 1103] 1104] 1105] and a full double-pole 
approximation )106] 1107] 1108] 1109] [60] are available in the literature. 
At the production threshold, additional large corrections can arise from the Coulomb 
singularity [110] llll] . For on-shell particles, the cross section factorizes according to 

a vr 

<7C = TT^ <^Born, 
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iVlcioo 


VV 1(1 ILL 


xt 


183.7 GeV 


0.077 GeV 


xt 


415.4 GeV 


3.1 GeV 



Table 3.1: Chargino masses and widths for the SUSY parameter set SPSla'. 



where P denotes the relative velocity of the produced particles. This expression diverges 
for /? — > 0. At threshold, the Coulomb singularity needs to be resummed using effective 
field theories. If the produced particles are off-shell, the singularity is cut off at a 
relative velocity (3 oc \JV jra. For the masses and widths of charginos given in Table 
13. H this leads to corrections < 1 % at threshold. Similarly, for W pair production 
\y\-1\ 1113] 11141 [TT5] and slepton production, |116j . corrections are in the percent regime. 
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4 Inclusion of NLO corrected matrix 
elements in WHIZARD (fixed order) 



4.1 Monte Carlo (event) generators 

In general, Monte Carlo techniques make use of random numbers to numerically deter- 
mine values of integrals or, given a probability distribution, simulate the outcome of 
physical events. For a general introduction of Monte Carlo techniques in particle and 
especially collider physics, see e.g. [in \ [57 1 1118] . 



4.1.1 Monte Carlo integration 

In Monte Carlo integration, the general idea is to use 

f{x)dx PS (X2 - xi) (/(x)), 

where {f{x)) is determined by the averaged value of N random calls of x: 

1 ^ 



i=l 



According to the central limit theorem for large numbers, the error is then oc 1/yN. 
The A^-dependent error can be decreased by importance sampling, where more values of 
X are chosen in regions where f{x) is largest, or similar techniques. 
For the numerical calculation of cross sections with n final particles, we need to integrate 



(4.1) 



where 



n 



(2 7r)^5(^) 



is the n-dimensional final state phase space. In Monte Carlo programs, the matrix 
element M is either coded manually or generated by some (internal or external) au- 
tomatic matrix element generator. Examples for external programs are CompHEP [119j . 
MadGraph |120j . or O'Mega [121j . The integral depends on 3n — 4 independent variables. 
Of course, the multi-dimensional integration of phase space is non-trivial and equally 
requires refined techniques. Differential cross sections can be obtained accordingly. 
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4.1.2 Event generation 

A physical event is defined by the specification of the n four-momenta of the final state 
particles, which require the generation of 3 n — 4 random numbers. In a straightforward 
Monte Carlo integration of a (14. ip . all events have the same a-priori probability. To 
obtain the final result a, they are weighted with the corresponding differential cross 
section da j (fl d-pi). 

A Monte Carlo Event generator, in contrast, should generate events according to their 
actual probability. This can be achieved by adapting the a-priori probability to the 
physical distribution or applying a hit-and-miss technique where each event is assigned 
the corresponding probability Pi = da [i-^/ da max which is compared with a random 
number r between and 1 and kept if Pi > r. Notice that this requires that dai > 0. 
Although this condition is always fulfilled in leading order, NLO calculations might cause 
problems for certain points of phase space [57] ; cf . Sections 14.51 and 15.11 
The events generated by the Monte Carlo program provide the same information as 
experimental data and can be analyzed using the respective detector simulation and 
analysis tools. Note that this equally allows for plotting of one- or multidimensional 
partial distributions, correlations, etc. without any further analytic calculation. 

4.1.3 WHIZARD 

WHIZARD [52j is a universal Monte Carlo event generator for multiparticle scattering 
processes. It interfaces several matrix event generators such as CompHEP , MadGraph , 
and O'Mega . In addition, it includes initial state radiation, beamstrahlung using the 
program CIRCE [122j, and fragmentation and hadronization routines from pythia [55j. 
It is designed as a 2 ^ n particle event generator. In one call, several processes can be 
combined such that background studies are simplified. Similarly, the results of several 
matrix element generators can be compared. For SUSY processes, this has recently 
been used for an extensive comparison [123j . Currently, it includes the Standard Model, 
MSSM, little Higgs models, and non-commutative geometry models. Similarly, it allows 
for user-modified spectra, structure functions, and cuts. 

4.2 Calculating NLO matrix elements using FeynArts and 
FormCalc 

FeynArts [1241 [86] and FormCalc |81j are Mathematica- and Form-based programs for 
(higher-order) matrix element generation and the calculation of the respective total and 
differential cross sections. It includes the SM, the MSSM, and can be extended to any 
model desired by the user. Furthermore, it can generate Feynman diagrams in a Latex 
or postscript format. Both programs use LoopTools |81j for the calculation of n-point 
functions and other loop-related quantities. We will quickly discuss both FeynArts and 
FormCalc and refer to the respective manuals [80[ [85] for more details. 
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FeynArts 



FeynArts is a purely Mathematica-based program. For a given number of in- and out- 
coming particles and loops, it first generates the corresponding general topologies. After 
choosing a physical model and specifying the in- and outgoing particles, all amplitudes 
for the specified process are generated and given analytically (depending on the process, 
the output might be quite complex). The user can then apply numerous specifications, 
e.g. diagram selections or renormalization conditions. FeynArts equally creates the 
Latex or postscript output for all created or specifically selected diagrams. 



FormCalc is a Form |125l I126j -based program with a Mathematica interface. It gener- 
ates a Fortran code corresponding to the FeynArts amplitudes. The resulting program 
numerically integrates the total or (angular) differential cross section for the correspond- 
ing process. Currently, it contains the kinematics for 1 — > 2, 2 ^ 2, and 2^3 particle 
reactions. In addition, it provides an easy input for e.g. mSugra parameters, loops over 
the cm energy or model parameters, or energy and angular cuts. It equally allows for 
different choices of multi-dimensional integration routines. The integrations are carried 
out according to the techniques described in Section 14.1.11 

Technically, the generated code consists of different process-dependent or independent 
modules. They contain e.g. routines for general features of the process, the process kine- 
matics, the code for the Born and one-loop matrix element and, for SUSY processes, code 
for the SUSY spectrum-generation. The compilation creates libraries for the calculation 
of the matrix element (squaredjne . a), the renormalization constants (renconst . a), and 
kinematics-dependent variables (util.a) which are linked to the main executable. Fur- 
thermore, the LoopTools library (libooptools . a) has to be included. 

4.3 Inclusion of the fixed order NLO contribution using a 
structure function 

In WHIZARD, there is an interface for arbitrary structure functions f{xi,X2) that can be 
convoluted with the cross section according to 



where Xj is the beam energy fraction. f{xi,X2) can be the sum of two uncorrelated 
structure functions (one for each incoming beam) or a correlated structure function for 
two incoming beams. For more details, cf. |52j . 




FormCalc 





) 
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We can therefore implement the fixed-order one-loop result atot (Eq. (j3.1ip ) using 
crtotis,ml) = j dx fesixi, X2;AEy,Aey,^) j dT2\McE{s,xi, X2;ml 



,2m2 



+ / d^3|^^2^3(s)|^ (4.2) 

with the structure function 

/efr(xi, X2] AE^, A9y, ^) = 6(1 - xi) 5(1 - X2) 

+ S{1 - xi) f{x2; AO^, ^) ^(xo - X2) 

+ /(xi; A0^, ^) 6{l - X2) 9{xo - xi) (4.3) 

and the effective squared amplitude 

\Mes{s,Xi, X2;ml)f = [l + /soft(AE^, A^) ^(xi, Xa))] \MBorn{s)\^ 

+ 2Re [MBoTn{s)Mi.ioopis, m^)] ^(xi, X2) (4.4) 

with 6'(xi,X2) = 0{xi — xq) 0{x2 — xq). The parameter 

2AE 

xo = 1 



separates the hard from the soft photon region. In addition to the introduction of the 
structure function /eg, we replace the Born matrix element as computed by the matrix- 
element generator, O'Mega, by the effective matrix element A^e// (14. 4p . The latter is 
computed by a call to the FormCalc-generated routine. 

Technically, (j4.2|) is implemented by splitting the structure function into four different 
regions in the (xi,X2) space: 



• Xi > Xo, X2 > Xq : (soft-soft) 

This corresponds to the region where /e// = (5(1 — xi)5(l — X2). Both photons 
are in the soft photon regime; we set Xj = 1 and calculate the Born+virtual+soft 
cross sections: 

CTBorn+Crvirt+CTsoft = j dT2 { (1 + /so/t) Born (s) | ^ + [2Re (A^Born(s)* A^l-loop(s))] } • 

All matrix element contributions in this region as well as the soft photon factor 
are generated by a call to FormCalc. As we are mapping two 5 functions to a finite 
(xi, X2) region, we need to divide the result by a normalization factor N such that 
4r dxi dx2 = 1. 
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• Xi > xq , Xj < xq (soft-hard) 

2 

This corresponds to the region where /g// = S{1 — Xi) /{xj; AOj,^). One of 
the photons is in the soft regime, while the other one is in the hard regime. For 
the soft photon, x is again set to 1. The hard-cohinear contribution of the hard 
photon is then generated according to 

^ /(/-W.S„J..») + 

= dxTix-Ae,,^) dT,\M'^,^^{xs,mir- 

Here, the matrix elements can be calculated using either FormCalc or 'Mega . The 
helicity-dependent form factors are implemented by a modification of the density 
matrix. In general, 

\M\' = MHX)pxx'M{X'), 

where pw is the helicity-dependent density matrix for the incoming particles. For 
the inclusion of it is modified such that 

P'w = f^^X,X' PX'X' + f~Sx-X' PX'X', 

where we assumed that A, A' can take the values of it 1 as in fermionic cases. Tak- 
ing the mapping of the 5 function into account, we again have to divide the result 
by a normalization factor. 

• xi < xq , X2 < xq (hard-hard) 

Both photons are in the hard regime. In principle, this region would describe the 
radiation of two hard photons. As this is a second order effect, we set /e// = in 
this region. 

As the soft photon approximation requires x ~ 1, we have to artificially split the 
interval [0 : 1] such that x > xq is reached sufficiently often. This can be done by the 
introduction of an additional x^d and the projection 

[0 : Xdei] — > [0 : xo] ; [xdei : 1] — > [xq : 1]. 

We then have to multiply all results by the Jacobian of the transformation. 

The hard, non-collinear part is added by a separate run of WHIZARD as given by eq. (j3.10p 

with the explicitly generated matrix element M2^3 applying the respective AE^, A6^ 
cuts. However, in WHIZARD it is equally possible to combine different processes in one 
run which then reproduces atot as given in Eq. (j3.1ip . 
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Implementing this algorithm in WHIZARD, we construct an unweighted event generator. 
With separate runs for the 2 — > 2 and 2^3 parts, the program first adapts the phase 
space sampling and calculates a precise estimate of the cross section. The built-in rou- 
tines apply event rejection based on the effective weight and thus generate unweighted 
event samples. 

On the technical side, for the actual implementation we have carefully checked that 
all physical parameters and, in particular, the definition of helicity states are correctly 
matched between the conventions [77] used by O'Mega /WHIZARD and those used by 
FormCalc (cf. e.g. |127j ). These differ by a complex phase in the definition of the two- 
component helicity eigenstates. 

4.4 Technicalities of the implementation 

The inclusion of the FormCalc generated NLO contributions described in Section 14.31 
implies the modification of the WHIZARD as well as the FormCalc code. Here, we just 
sketch the general approach and refer to [128j for more details. 

We used the FormCalc code for chargino production at NLO available from [64^ I84j. 
WHIZARD already includes calls to other external programs, such as CompHep, Mad- 
Graph, and O'Mega. Here, libraries are created from the respective programs, which are 
then linked to the WHIZARD executable. The same is done here for the inclusion of the 
NLO FormCalc matrix element. 

We therefore modify, in the FormCalc code, those modules calculating the matrix ele- 
ment as well as the soft photon factor (to create a WHIZARD -compatible in- and output) 
and the kinematics (in FormCalc, the routines setting the kinematics equally set the 
soft photon factor and give values to internal variables needed for the matrix element 
evaluation). In addition, we create a new interface setting the SM and MSSM variables 
to the values used by WHIZARD . We can therefore make use of the internal WHIZARD data 
reading routines and can equally use input in the Les Houches accord format [129j . A 
library libf ormcalc . a is then created and linked. The library equally contains the 
(modified) FormCalc libraries for the squared matrix element (squared_me . a), some 
kinematic variables (util.a), and the renormalization constants (renconst . a). 
In WHIZARD, we add a routine set_pars_f ormcalc setting the SM and MSSM values used 
by WHIZARD in the FormCalc routines as well as a call to the (NLO)-generated matrix 
element subroutines SquaredME which depends on the cm energy and the four-vectors 
of the outgoing particles. 

To summarize, we 

• modify the file 2to2.F from FormCalc, containing the 2 — > 2 kinematics, for in- 
and output of the WHIZARD generated particle momenta, 

• modify the file squared_me . F such that the helicity-dependent subroutine 
SquaredME (result , spins, reset, matrix_Born, & matrix_loop, old_angles, 
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old_s) can explicitly be called from WHIZARD, 

• generate a new file transer . f for SM and MSSM value transfer, 

• create the library libf ormcalc . a from the modified code, including also the un- 
modified libraries renconst . a and util . a, 

• link libf ormcalc . a and libooptools . a to WHIZARD. 
In WHIZARD , we 

• add the respective routine set_pars_f ormcalc for variable transfer in the process- 
file, 

• call the subroutines setting the case-to-case kinematic variables SetEnergy and 
SquaredME each time we need the FormCalc one loop matrix element and soft 
photon factor. 

In addition, we included checks guaranteeing that the FormCalc loop-related routines 
are only called when necessary. For example, ^/s dependent quantities are recalculated 
only if actually changes. 

The difference between the helicity bases for fermions used in WHIZARD and FormCalc (cf. 
Section 14. 3p needs to be accounted for by the introduction of an additional phase of the 
matrix element. For the Born contribution, this has been checked for all possible helicity 
combinations. 

The soft and collinear cuts A E^, A6^ are added as variables in the file cutpars . dat in 
the results subdirectory of WHIZARD. 

4.5 Drawback of the fixed-order method 

For any fixed helicity combination and chargino scattering angle the differential cross 
section is positive if we include all contributions defined in Eq. (jS.lip . If the integration 
and simulation is split into a 2 — > 2 and 2^3 part as implied by atot P-Sp . however, 
the fixed-order approach runs into the well-known problem of negative event weights \88\ 
\130\ [57] : The effective 2 — > 2 matrix element in the soft-soft integration region, 

\Meff\^ = {l + fsoft)\MBorn{s)\^ + 2Re{MBorn{s)MU,,p{s)), (4.5) 

is no longer positive definite if AE^ becomes sufficiently small. If we lower the cutoff, 
this becomes negative within some range of scattering angle. We will investigate this 
effect closer and concentrate on the dominant helicity contribution with 

Ae+ = 1 , Ae- = -1 , A^- = 1 , = -1, (4.6) 

and consider the A, AE^ dependencies of |A^e//P- 
In the effective squared matrix element (14. 5p . we have 

fsoft = fsoft{X,^E^) ; Ml-loop = Ml-loopW- 
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The AE"^ dependence enters only in the soft factor 

fsoft tx In 

In the combination of virtual and soft photons in Mf,ff, the log A dependence cancels 
exactly. 

As an example, we consider the angular behavior for a fixed helicity state with = 
0.005 (high cut) and —7^ = 0.0005 (low cut), respectively. Figures BTT] and show all 

contributions to |A^e//P in these cases. M-Bom and A^i_;oop(A) are AE^ independent 
and therefore do not change (for a fixed value for A). With the high cut, we obtain 
IMeffl"^ > for all values of 9, while for the lower cut, |A^e//P < for cosO < 
0.4. In these regions of phase space, is small enough such that the virtual photon 
contributions are not sufficiently canceled by soft real photon contributions. Figure 
14.31 shows the behavior of the |A^e//P for the helicity combination (|4.6|) as well as the 
subdominant contribution with A;^- = A^+ = 1. 

If we insist on a positive weight Monte Carlo generator, an ad-hoc solution for the 
fixed 0{a) contribution would now be to set -A^e// = in the respective regions of 
phase space. However, for too low soft cuts, this leads to wrong results for the total and 
differential cross section, cf. Figure [131 An alternative approach uses subtractions in the 
integrand to eliminate the singularities before integration |13H 11321 11331 \T3^ 1101^1135] . 
The subtracted pieces are integrated analytically and added back or canceled against 
each other where possible. However, the subtracted integrands do not necessarily satisfy 
positivity conditions either. 

Alternatively, we can include negative event weights. Such event samples are not a 
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Figure 4.2: |A^e//P contributions for ^ = 5 • 10-^ LO: Born, NLO: |A^e//P, v 
denotes the virtual and s the soft contribution 




Figure 4.3: Effective squared matrix element (arbitrary units) for €~ — > x{yi[ as a 
function of the polar scattering angle Q at = 1 TeV. Left figure: Helicity 
combination — h H — ; right figure: — V ++. Solid line: Born term; dashed: 
including virtual and soft contributions for A£^^ = 10 GeV; dash-dots: same 
with A^;^ = 0.5 GeV. 



37 



160 



150 



140 - 



130 



fix 



= 1 TeV 



LO 



0.1 



10 



100 



AE-, [GeV] 



Figure 4.4: Total cross section as a function of the energy cutoff AE^ using different cal- 
culational methods: (green, dashed) = fixed-order semianalytic result using 
FeynArts/FormCalc; (red, dotted) = fixed-order Monte-Carlo result us- 
ing WHIZARD; semianalytic result for lowest numerically reachable AE^. Ef- 



fects from setting |A^e//| 



where it becomes negative are visible for 



AE^I ^/s < 10 . Statistical Monte-Carlo integration errors are shown. For 



the Monte-Carlo results, the collinear cutoff has been fixed to A^^ 



1°. 



The rise of the fixed order result reflects the breakdown of the soft photon 
approximation; see Section [4.61 
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Fi gure 4.5i Effective squared, matrix element cind. squared. NLO contribution |yVf^_|_5| as 
a function of the polar scattering angle at ^/s = 1 TeV. Helicity combi- 
nation h H — ; Solid line: Born term; dashed: including virtual and soft 

contributions for AE^ = 5 GeV; long dashed: same with AE^ = 0.5 GeV; 
dash-dotted: A = 5 • IQ-^ GeV 



possible outcome of a physical experiment and imply further modification of the detector 
and analysis tools. 

In general, experimental resolution at the ILC can well reach = 10^"^ and lower 

values. From Figure 14. 5t we see that, for a fixed A, the 1-loop contribution |.A^]^_;Qop| 
is nearly the same order of magnitude as |A^e//P- However, the soft contribution here 
is treated with the infrared cutoff A; setting A ^ will always lead to the well-known 
infrared divergence (Tso/t oo. In principle, this can be canceled by the terms describing 
the emission of two soft photons equally regulated with the photon mass A. Therefore, 
if we want to construct an event generator reaching the experimental cut requirements, 
we should take second and higher-order contributions into account. 
For low enough soft energy cuts, the problem of negative event weights and negative 
soft cross sections remains as long as only finite order photon emissions are taken into 
account ([57J and references therein). This signals a breakdown of perturbation theory 
in this region of phase space and requires the inclusion of summed ISR contributions 
discussed in Section 13. 7[ We will present a method to combine this with the fixed-order 
corrections amrt in Chapter [3 
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4.6 Results 



4.6.1 Cut dependencies 

In the kinematical ranges below the soft and colhnear cutoffs, several approximations are 
made. In particular, the method neglects contributions proportional to positive powers 
of AE'^ and A6^, so the cutoffs must not be increased into the region where these effects 
could become important. On the other hand, when decreasing cutoffs too much we 
can enter a region where the limited machine precision induces numerical instabilities. 
Therefore, we have to check the dependence of the total cross section as calculated by 
adding all pieces and identify parameter ranges for AE^ and A6^ where the result is 
stable but does not depend significantly on the cutoff values. 
In the following, we will compare 

1. fixed order: The implementation of atot according to Eq. ()4.2p in WHIZARD with 
both soft (A£^^) and collinear (A^^) cuts, 

2. semianalytic: the NLO calculation presented in [63] using FormCalc in combi- 
nation with the 2 — > 3 part from WHIZARD . Here, only a soft cut AE^ is applied. 

Both programs use the same routine for the calculation of Mvirt and fsoft as well as the 
same SM and MSSM input parameters. Therefore, differences are due to the use of the 
collinear approximation and implementation differences described in Section [4.11 
Throughout this section, we set the process energy to y/s = 1 TeV and refer to the 
SUSY parameter point SPSla'. All 2 — > 2 and 2 — > 3 contributions are included, so the 
results would be cutoff-independent if there were no approximations involved. 

Energy cutoff dependence 

In Figures [4.41 and WM. we compare the numerical results obtained using the semianalytic 
calculation with our Monte-Carlo integration in the fixed-order scheme. The semiana- 
lytic result is not exactly cutoff-independent. Instead, it exhibits a slight rise of the 
calculated cross section with increasing cutoff; for AE^ = 1 GeV (10 GeV) the shift 
is about 2 %o (5 %o) of the total cross section, respectively. This is an effect of the soft 
photon approximation, where the energy fraction x of the incoming electron/ positron 
is set to 1 in the soft regime. Therefore, for x ~ 1 the error is 0{AE^/^/s) with respect 
to 

^Borni cf. Figure 14. 6[ For AE^/^/s < 10 ^, we run into numerical problems with 
the exact 2^3 process. Otherwise, the errors of the semianalytic calculation are in 
the per mille regime and smaller. 

The fixed-order Monte-Carlo result agrees with the semianalytic result, as it should be 
the long as the cutoff is greater than a few GeV. For smaller cutoff values the 

Monte-Carlo result drastically departs from the semianalytic one because the virtual 
correction exceeds the LO term there, and therefore the 2 — > 2 effective squared matrix 
element becomes negative in part of phase space as discussed in Section 14.51 There, the 
integrand is set to zero. 
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Figure 4.6: Total cross section as a function of the energy cutoff Ai?^ using different 
calculational methods: fixed-order semianalytic result; fixed-order Monte- 
Carlo result using WHIZARD. Relative change with respect to atot for ^E^^ = 
0.1 GeV, no angular cut 



Collinear cutoff dependence 

The collinear cutoff separates the region where, in the collinear approximation, 
higher-order radiation is resummed from the region where only a single photon is in- 
cluded, but treated with exact kinematics. We show the dependence of the result on 
this cutoff in Fig. 14.71 We see that, for A ^ < 1°, the collinear approximation holds. For 
A^ = 1°, its errors are less than per mille. For larger A0, the collinear approximation 
breaks down. Similar results are found in [84^ I66j. 

Photon mass dependence 

The infinitesimal photon mass A is used in the FormCalc matrix element calculation for 
the regularization of infrared divergencies. The effective matrix element |A^e//P (Eq. 
14. Sp should then be independent of A. Numerically, this has been tested for 

5 X 10"^ GeV < A < 10^° GeV 

for the FormCalc integration routines. In these regions, while the photon mass remains 
a parameter in the matrix element code, the result does not numerically depend on it, 
regardless which method has been used. It is of course more meaningful to choose small 
photon masses. 
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Figure 4.7: Total cross section dependence on the collinear cutoff A6^ using different 
calculational methods: sa (red, dotted) = fixed-order semianalytic result fix 
(green, dashed) = fixed-order Monte-Carlo result atot (|3.1ip : The soft cutoff 
has been fixed to AE^ = 10 GeV, ^/s = 1 TeV 



4.6.2 Total cross section 

Fixing the cutoffs to AE^/^ = 5 • 10"^ AO = 1°, we can use the integration part 
of the Monte-Carlo generator to evaluate the total cross section at NLO for various 
energies. They exactly reproduce the semianalytic results modulo the A i?^-dependent 
%o errors discussed in Section 14.6.11 For a discussion on the general behaviour of the 
fixed NLO cross section, cf. Chapter [3l The numerical results qualitatively agree with 
[6^ [66] . However, we did not compare all SUSY parameters used here. 

4.6.3 Event generation 

The strength of the Monte-Carlo method lies not in the ability to calculate total cross 
sections, but to simulate physical event samples. We have used the WHIZARD event 
generator augmented by the effective matrix element (14. 5p and structure function ()4.3p 
to generate unweighted event samples for chargino production. 

To evaluate the importance of the NLO improvement, in Figures 14.81 and 14.91 we show 
the binned distribution of the chargino production angle as obtained from a sample of 
unweighted events corresponding to 1 ab~^ of integrated luminosity for ^/s = 1 TeV 
and the SUSY parameter point SPSla'. With cutoffs AO-^ = 1° and AE^ = 3 CeV we 
are not far from the expected experimental resolution, while for the fixed-order approach 
negative event weights do not yet pose a problem. 

The histograms illustrate the fact that NLO corrections in chargino production are not 
just detectable, but rather important for an accurate prediction, given the high ILC 
luminosity. The correction cannot be approximated by a constant proportionality factor 
K between the leading and next-to-leading order cross section (ET-factor) but takes a 
different shape than the LO distribution. The correction is positive in the forward and 
backward directions, but negative in the central region. For comparison. Figure also 



42 



2x10^ 





1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 




e"e+ XiXi 






evt/bin 


— -. 

j 1 

1 i - 


- 


= 1 TeV 


: i 1 

r-l" \ - 


- 


JC = 1 ab"^ 


^ 

i ' 

r " "' 

1 


f 


■ ■■) 1 1 1 1 


■ i 

1 — 

1 1 1 1 1 1 1 1 1 1 1 1 1 



-1 -0.5 0.5 1 

cos 9 

Figure 4.8: Polar scattering angle distribution for an integrated luminosity of 1 ab~^ at 
^/s = 1 TeV. Total number of events per bin, Born (dotted, black) and fixed 
order (dashed, green) result 

shows the 1 a error from the Born result. 
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5 Resumming photons 



5.1 NLO and resummation at LEP and LHC 

The principle of analytic photon resummation has already been discussed in Section [3.71 
We will here shortly sketch the implementation of photon resummation in Monte Carlo 
generators at LEP as well as NLO matrix elements for hadronic machines and refer to 
the literature [571 UlSj for more details. 

Photon resummation at LEP 

The problem of negative event weights discussed in Section H3] was. at LEP times, known 
as the fco problem. When slicing the Monte Carlo generator phase space in a 2 — > 2 
(Born, virtual, soft) and a 2 — > 3 (hard) part as in the fixed order approach, event 
weights in the 2 — > 2 part can become negative for ko 0. This holds true in any 
finite-order treatment of soft real and virtual photons. The solution is to sum photons 
up to an infinite order. In ^7\, different methods are proposed in order to include 
exponentiation in Monte Carlo event generators: ad hoc-exponentiation which estimates 
the energy loss due to photon radiation and modifies the Monte Carlo code accordingly, 
the use of (positive definite) structure functions convoluted with the Born cross section 
at a reduced cm energy, and YES exponentiation [136j which correctly simulates the 
hard part, while using a structure function for soft and virtual photons. In Section \37l\ 
and in the following, we will use the structure-function method proposed by Jadach and 
Skrzypek [75]. Note, however, that none of these approaches includes non-photonic NLO 
contributions as amrt (Eq. 13. 2p . 

NLO and parton showers for LHC 

Recently, several computer codes attempt to include exact NLO matrix elements in 
Monte Carlo generators containing partonic showers, which are the QCD equivalent of 
the infinitely summed up photon contributions described above. These are basically 
directed to experiments at hadronic machines such as the LHC, cf. [TU (721 11181 173] . A 
similar code treating jet production at e"*" e~ colliders was proposed in [69l[70]. Basically, 
the NLO contributions are added by hand, while the corresponding first order parts from 
the parton showers are subtracted. In the following, we will pursue a similar approach 
for the inclusion of electroweak NLO corrections. 
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5.2 Resummation method 



The shortcomings of the fixed-order approach described in Chapter H] are associated 
with the soft-colhnear region i?^ < AE^, 0^ < A6^, where the appearance of double 

logarithms ^ In ^ In 6^ invalidates the perturbative series. However, in that region 
higher-order radiation can be resummed [94^ 11371 I138j . The exponentiated structure 
function /isr ()F.14p that resums initial-state radiation, 

o-Born+iSR(s, A6'^,mg) = y dx fisRix; AO^ , ^) j dT2\MBQrn{xs)\'^ , 

includes photon radiation to all order in the soft regime at leading-logarithmic approxi- 
mation and, simultaneously, correctly describes collinear radiation of up to three photons 
in the hard regime. It does not account for the helicity-flip part /~ (|F.5|) of the fixed- 
order structure function; this may either be added separately or just be dropped since 
it is subleading. 

We now combine the ISR-resummed LO result with the additional NLO contributions 
f^virt given by Eq. (j3.2p . To achieve this and avoid double counting, we first subtract 
from the effective squared matrix element, for each incoming particle, the contribution 
of one soft real and virtual collinear photon that is contained in the ISR structure func- 
tion. The soft photon has already been accounted for in the soft-photon factor, while 
the virtual part is contained in the interference term. Then 

|7Wefr(s;A^^,Ae^,m2)|2 = + _ 2/,,ft,isR(A^^, A0^, ^)] |7WBorn(s)P 

+ 2Re [MboUs) -Mlioop(^> m2)] (5.1) 

with s being the cm. energy after radiation and /soft,iSR the integrated 0[a) contribution 
of fiSR- This expression contains the Born term, the virtual and soft-collinear contri- 
bution with the leading-logarithmic part of virtual photons and soft-collinear emission 
removed, and soft non-collinear radiation of one photon; it still depends on the cutoff 
AE^. Convoluting this with the resummed ISR structure function. 



o-v+s,isr(s, AS^, AO^, m, 

= j dxifisR{xi;Ae^,^) j dx2 /isr(x2; A0^, ^) j dTs |7Weff(s; A^^, A^^, mi 



2m2 



(5.2) 

we obtain a modified 2 — > 2 part of the total cross section. 

In this description of the collinear region, there is no explicit cutoff AE^ involved, 
and collinear virtual photons connected to at least one incoming particle are included. 
The cancellation of infrared singularities between virtual and real corrections is built-in 
for collinear photons. The main source of negative event weights is eliminated, and we 
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obtain a better behaviour for the integrand such that smaller energy cuts can be ap- 
plied; cf. Section [5.4. 11 Furthermore, the method is exact in 0{a) leading log terms. In 
the following, we will consider the description of one and more photons resulting from 
cJv+sjsR in more detail. 

In the resummation method, the emission of real and virtual photons is described in 
various ways: the soft approximation (cf. Sec. 13. 3p . initial state radiation (cf. Sec. 13. 7p . 
virtual contribution from interference term (cf. Section 13. 2|) , and real emission given by 
exact (hard non-collinear) matrix element M2—>3 (cf. Sec 13. 5p . We will now consider 
the respective description of photons resulting from the resummation method in different 
points of phase space. Here, we go to infinite order for the photon emission from one 
incoming particle and to second order for the simultaneous photon emission from two 
incoming particles. 



Radiation off one incoming particle 

We first consider photon radiation from one particle only and use the exponentiated 
electron structure function fisR (|F.14p . In the following, we ignore the emission of the 
second photon contained in fsoft ^-nd only subtract the contribution of fisR for one 
photon. The corresponding effective matrix element for the radiation of only one photon 
is 

|7wS(s; AE^, A9^, ml)\^ = [l + /soft(^) - /soft,iSR(Ai?^, A9^, ^)] |A^Born(s)|' 

+ 2Re [MBornis) A^lioop(^> A^ ml)] (5.3) 

and the total cross section 

4+s,iSR(^.^^7.A^7'"ie) = / dxfisRix;A9^,^) J dr2\M'^J^{x,s;AE-y,A9^,ml)\\ 

(5.4) 

According to the exponentiation principle discussed in Appendix IF. 3. 21 we can decom- 
pose (J) fisR in factors F^^\Fi such that 
»i n 

/ D^^{x)dx = ^ fW(xo) + 0(7?"+^), 

n 

D^^{x) = ^Fi(x) + 0(??"+i) x<xo, (5.5) 
/=o 

where F^'^\xq) is the Oirf") contribution to the integrated exact solution of the evolution 
equation ()3.13p in the soft limit and Fn{x) is the exact 0{r]'^) perturbative solution to 
it. In the structure function fisR (|F.14p used in this work, F„ has been calculated up 
to n = 3. Using the scale = {A9pq)'^, where po is the energy of the electron, these 
factors exclusively describe collinear photons. We have 

^(0) ^ 1 ^ ^(1) ^ f^^f^^jsR ; i^o = 0. (5.6) 
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|A1e//P mixes different orders of r]. For a fixed order n, the photon contributions to 
4?s,iSR are given by 



where 



= F^''\xo)aBorn(s) + / Fn{x) aBornix, s) dx + F^""-^^ a{s) 

Jo 

+ / Fn-i{x)a{x,s)dx-F^'\xo)F^''-'\xo)aBorn{s) 
Jo 

-F^^\xq) I Fn-l{x)aBorn{x,s)dx, (5.7) 
^0 

a = y" dr2fsoh\MBoTn{x,s)f + 2I{B[MBornix,s)Ml.loop{x,s)] 

contains all contributions from virtual photon emission according to the exact A^i-ioop 
calculation as well soft photon emission in the soft approximation. 

We will now consider the results for soft-collinear photons, distinguishing between multi- 
ple photon emissions where the last emitted photon does or does not obey the transverse 
momentum ordering k±^n-i ^ ^±,n- 

• A;_L ordering obeyed 

If the first n — 1 emitted photons are soft (i.e. their total energy is lower than 
AEy), this part of is given by 

F^^\xo)aBorn{s) + F(--'^ais)-F('\xo)F(--'\xo)aBorn{s) = F(--'^a{s), 
and by 

/ F„(x)-F«(a;o)F„_i(x) ^ Born{x ^ S) dx \ Fjx—\{x)o{x^S)dx 



Jo 



Fn-i{x) a{x, s)dx 



otherwise. The last soft-collinear photon is described by the soft photon approxi- 
mation, all others by fisR- 



• no k± ordering 

The terms oc F^^\ F„, which describe the emission of n A;_L-ordered photons, are 
missing in This results in 

(d{s) - F^^\xo)aBorn{s)) + ^ i^n-l(x) (d{x,s) - F^^\xo) a Born{x , s)) dx. 

(5.8) 
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Figure 5.1: Example for mixed real and virtual photon emissions described by /isr- All 
photons are collinear and fcj^-ordered 

We only obtain differences between the exact and the leading log contribution for 
the last radiated photon, multiplied with the fisR contribution for the first n — 1 
photons. 



If only hard-collinear photons are emitted, they are all described by fisR- For virtual 
photons, the same relations hold as for soft real photons: in the case of kx_ ordering of 
the last photon, it is given by the contribution to Mi_ioop- Otherwise, the emission of 
the virtual last photon is again described by difference terms according to Eq. (jS.Sp . 
Note that the radiation of a non-collinear soft photon as well as a non-coUinear virtual 
photon is not touched by the subtraction mechanism. Here, the respective contribution 
to fj^^^ is simply 



JO 

where Orest contains all non-collinear virtual and non-photonic contributions. 
In general, fisR describes a combination of real and virtual photon emissions, e.g. graphs 
given in Figure EUl as the electron splitting function Pge (|3.12|) combines the radiation 
of both a real and a virtual photon in order to cancel the IR divergence. Both soft real 
and virtual photons are contained in the soft contributions to the factors F^^^ and -F„. 
For example, we can split F^^^ according to 



This describes the emission of n collinear photons in an arbitrary combination of soft 
and virtual contributions. The only requirement is the ordering of their transverse 
momenta and X^^eal 7 -^7 — ^-^7- 

Radiation off two incoming particles 

Figure 15.21 shows the diagrams contributing to the emission of two photons off two in- 
coming particles for the s channel 7, Z exchange (we omitted diagrams resulting from 
the crossing symmetry for the radiated photons). As discussed in the last section, virtual 
photons are contained in the soft parts of jisR (cf. Eq. ()5.10p ). We now consider the 




(5.9) 



f(")(xo) = (F(i))- = (F,„^, + F,,rt)" 



(5.10) 
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Figure 5.2: Diagrams contributing to the radiation of two photons, either from one in- 
coming particle, or both incoming particles, for the s channel 7, Z exchange. 
Missing are diagrams from crossing symmetry of radiating photons. The t 
channel iJ-exchange behaves analogously 

contributions to (Jmrt+ISR (15. 2p from up to two photon emissions. In the remainder of 
this section, a = a Born- 

In general, the photonic 0(a") contributions to o"v+s,isr are given by 

V / dxi dx2Fi{xi)Fn-i{x2)cr{xiX2s) + / dxiF,(xi)F("~')(x2)a(xis) 



+ / dx2Fn-i{x2)F^'\xi)a{x2s) + f("-')(x2) 



n' = n— 1 



" — " ^ / r-xo r-xo . 

+ V] ( / dxi / dx2Fi(x\)Fn'-i{x2) \a{x\X2s) -2F^^^ a{x\X2s) 
1=0 Vio Jo ^ 

+ " fixiFi(xi)F("'-')(x2) (<t(xis) -2FWa(xis 
+ ^ ° dx2Fn'-iix2) F^^\xi) (a{x2s)-2F^^^ a{x2S 
+ F«(xi)F("'^')(x2) (<t(s) -2F«ct(s))) . 



For n = 0, we obtain (recall Eq (|5.6|) ) the Born contribution (|3.1|) : 

xo rXQ nxo 

dxi / dx2 Fo{xi) Fo{x2) (t{xi X2 s) + / (ixiFo(xi) F*^°^(x2) (t(xi s) 
^0 Jo 

+ r (ix2Fo(x2)F(o)(xi)a(x2s) + F^^\xi) F^^\x2) a{s) 

Jo 
= ais). 
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Setting n = 1 gives the additional terms 



rxo rXQ rxo 

/ dxi / dx2Fi{xi)Fo{x2)a{xiX2s) + / dxiFi{xi) F^^> {X2) a{xi 
Jo Jo Jo 

+ ^ " dxiFo{xi)F^^\x2)a{xis) + F^^^xi) F^^\x2) a{s)^ + (xi ^ X2) 

rxo rxo , ^ 

+ J dxi J dx2 Fo{xi) Fo{x2) \^a{xi X2 s) - 2 F^^' a{xi X2 s) ^ 

+ ^ " dxiFo(xi) f(0) (x2) (<7(xi s) - 2 F^^) a(xi 
+ y dx2Fo(x2)F(o)(xi) (<t(x2s) -2F«CJ(X2S)^ 

+ (xi) (x2) (<t(s) - 2 a{si 



rxo rxo 

= / (ixiFi(xi) cr(xi s) + / (ix2-Fi(x2) ct(x2 s) + 
Jo JO 

where we used Eq. (j5.6p . the fact that F^^) only depends on xq such that F^^\x) = 
F'^^^xo), and 

-2FW = -FW(xi) - fW(x2). (5.11) 
for the subtraction terms in |A4e//P- 

We see that the n = 1 term corresponds exactly to the 2^2 part of the leading 
log 0{a) contribution of the fixed order calculation ()3.1ip . 

For n = 2, we obtain the additional terms 

(^F(2)(x2)fT(s) + ^ ° (ix2F2(x2)a(x2s)^ + (xi ^ X2j 

rxo r^o r-co 

+ / dxi / (ix2 -Fi(xi) Fi(x2) cr(xi X2 s) + / dxiFi(xi) F*-^-* (X2) (t(xi s) 
Jo Jo Jo 

dx2Fi(x2)FW(xi)a(x2s) + F«(xi) fW(x2) a(s) 

y dX2Fx (X2) (^(7(X2 S) - 2 F(1) fT(x2 s) j 

+ fW(x2) {a{s) - 2FW a(s))) + (xi ^ X2) , (5.12) 



where we again used Eq. (j5.6p . To really understand the Feynman diagrams to which 
the different contributions in ()5.12p correspond and the cancellations, we have to take 
Eq. (|5.1ip into account. We first consider the case of two photons radiating off the same 
incoming particle, e.g. the one depending on xi. The relevant contributions are given 
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by 

F(2)(xi)a(s) + / ' dxiF2ixi)a{xis) 



+ 



!•() 



(5.13) 



+ 



A comparison shows that this exactly corresponds to the n = 2 case of only 1 particle 
radiating off photons (cf. Eq. (|5.7p ). Therefore, the description of photon radiation 
given for this case also applies here. 

The mixed case is more complicated. The relevant terms are given by 

/ dxi / (ix2 Fi(xi) Fi(x2) cr(xi X2 s) + / dxiFi{xi) F^^\x2) cr{xi s) 
Jo Jo Jo 

f-xo 

+ / dx2Fi(a;2)F«(xi)a(x2s) + F«(x2) ct(s) 

dx2Fi{x2) ((7(2:2 s) -FW(xi)a(x2s)) + ^^^^(^2) (<t(s) - F«(xi) cj(s)) 

+ y ° dxiFi(xi) ((t(xis) -F«(x2)a(xis)) + fW(xi) ((t(s) - fW(x2) a(s)) . 

(5.14) 

If at least one of the radiated photons is hard, the relevant terms are given by (for e.g. 
xi hard and X2 hard or soft): 

/ dxi / dx2Fi{xi)Fi{x2)<j{xiX2s) + I dxiFi{xi) F^^\x2) (t{xi s) 
Jo Jo Jo 

+ j " dxiFi{xi) (a(xis) -FW(2;2)ct(xis)) 

f-xo f-xo f-xo 

= / dxi / dx2 Fi{xi) Fi{x2) cr{xi X2 s) + / dxiFi{xi) a{xi s), 
Jo Jo Jo 

i.e. we obtain the usual description (hard photon: leading log, soft photon: soft approx- 
imation, virtual photon: one- loop description from interference term). However, if both 
photons are soft, we have to closer investigate the phase space slicing. The soft-soft 
terms are given by 

F^'^\x2)a{s) + F^'^^xi)a{s) - F^^\xi) f'^^\x2) a{s). (5.15) 
We now split a into 

a = f(^) a +A F^^'^a + dresu (5.16) 

where o-rest contains all non-collinear virtual and non-photonic (ie, weak and SUSY) 
contributions, and AF^^^ is the difference between the soft approximation fsoft (|3.6p 
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and the first order contributions of the integrated leading log structure function fsoft,lSR 
(fRT8]l : 

A^^'^ = fl2ft - fsoftJSR. (5.17) 
In an ideal case, we would have 

(fW + AFW)(2:i)(F(i) +AF(i))(x2)ct. 

Up to O ((AF^^))^) terms, this corresponds to eq. ()5.15p . In this accuracy, both pho- 
tons are described by fsoft- The same holds for virtual-soft and virtual- virtual emissions, 
where the virtual part is (up to similar corrections) described by the interference term. 
Finally, we have to consider whether these contributions correspond to the photon- 
photon diagrams as given in Figure 15.21 We obtain 6 diagrams (for each diagram, 
there is one with the outgoing photons crossed). However, as we have two indistin- 
guishable photons in the final state, we obtain an additional factor ^ so that the ap- 
proximation above directly corresponds to the contributions we expect for the process 
e"*" e~ — > XX 77- The same holds of course for any convoluting of the Born cross 
section with J dxi f dx2 ^ [xi) ^ {x2) ■ 

As in the case of only one particle emitting photons (Eq. ()5.9I) ). (15. 2p also includes 
collinear photonic corrections to the Born/one-loop interference arest in leading log ac- 
curacy. The corresponding convolution with fisR is completely unaffected by the sub- 
traction mechanism, which only accounts for the soft parts. 

The complete result is supplemented by the 2 — > 3 part, 

'7tot,iSR(s,me) = fJv+s.isR + / dT^\M2^z{s)\'^ . (5.18) 

J AE-yAS-l 

A final improvement is to also convolute the 2 — > 3 part with the ISR structure function 
which defines 

a-tot,ISR+('S,?^e) 

= j dxifisR{xi;Ae.f,^) j dx2fisRix2;Ae^,^) 

X ( / dT2\M,s{s;AE^,Ae^,ml)\^+ [ dTs\M2^3{s)n ■ (5-19) 

This introduces another set of higher-order corrections, namely those where after an 
arbitrary number of collinear photons, one hard non-collinear photon is emitted. This 
additional resummation does not double-count. It catches logarithmic higher-order con- 
tributions where ordering in transverse momentum can be applied. Other, logarithmi- 
cally subleading contributions are missed. In 0{a), atot (|3.1ip is exactly reproduced 
by <Jtot,lSR'i only the helicity flip-part of the hard-collinear radiation is not taken into 
account (which can become important for completely polarized initial particle states). 
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For n = 2, crtotjSR ^ow also contains diagrams where the last radiated photon is hard, 
non-collinear, and obeying or not obeying the strong k± ordering. 

Leading and higher-orders: summary 

We can therefore summarize that crtot,iSR+ (15.190 reproduces atot (|3.11|) up to 0{a). 
For terms 0{a^) and higher, crtot,lSR+ contains all contributions of arest convoluted 
with fisR. for both incoming particles, i.e. the radiation of infinitely summed up soft 
or virtual collinear photons and up to 3 hard-collinear photons off all non-photonic and 
non-collinear virtual contributions of the interference term. The same holds if the last 
radiated photon is non-collinear and hard, as the contribution is then given by the 2^3 
part of (I5.19p . For the emission of a last soft or virtual collinear photon, we carefully 
have to check the contributions resulting from the subtraction. For 0{a'^), at least one of 
the photons is always described by the ISR structure function. But when the Born term 
is convoluted with the ISR function, there are also two-photon contributions described 
solely by the ISR. We have to distinguish between the cases where (i) the two photons 
are attached to the same or (ii) to different incoming particles. 
In case (i), we consider the three terms (cf. Eq. (|5.13|) ) 

Oia\sR - 0(a)isRO(a)fe| + 0(a)isR 0(a),oft. (5.20) 

The first term contains all pairs of collinear photons from the ISR, A;_L-ordered; the last 
term contains a first photon from ISR and a second one from the soft-photon factor or the 
interference term. The term in the middle is the subtraction to avoid double-counting of 
soft photons. Here both photons are from the ISR, the first one with arbitrary energy, 
the second one real soft or virtual. 

If the second of the considered photons is real soft or virtual, and both are /c_L-ordered, 
then there is an exact cancellation between the first two terms. For non /c_|_-ordered 
photons, the first term gives no contribution, and there is a cancellation between the 
second and third term, which results in a difference between the soft approximation/ 
interference term contribution expression and the ISR leading logarithmic approximation 
term given by AF(i) dSUD. 

In the case (ii), we write the terms schematically as (cf. Eq. (I5.14|) ) 

0(a)i,iSR0(a)2,iSR+0(«)i,iSR (o(«)2,.oft - 0(a)|^|R) + (o(a)i,soft - 0(«)i?isr) 0(«)2, 

Since there are always two different structure functions involved, /!;_|_ -ordering is absent, 
and after a cancellation of soft terms one is left with 

AF(i)(xi)0(a)2,isR + 0(a)i,isRAF(^)(x2) + 0(a) l,ISRC(a)2,ISR, 

which is up to the missing terms AF^^\xi) AF^^\x2) equivalent to a soft approxima- 
tion/ interference term description for both legs. 
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Remark: Soft approximation for matrix elements 

In the soft approximation as well as the exponentiation in structure functions, it is 
assumed that 

M{x > xo) ~ M{x = 1) , a{x > xq) ~ a{x = 1) (5.21) 
and therefore (cf. Eq. (|5.5|) ) 

/' dx fisdx) a{x, s) ^ ( [' dx fisnix)) a{s) = (1 + F^) a{s) + 0{rf) (5.22) 

J XQ \J XQ J 

for the exponentiation; fsoft is derived similarly. This holds true up to errors propor- 
tional to J^^ dx [da / dx)x = i fiSRi^) (1 — x). Some of these contributions, however, are 
accounted for in the actual WHIZARD implementation of the resummation method; cf. 
Section 15. 3[ 

To take them into account, we have to substitute /<jo/tO'(s') by 

fsoftcr{s') + i fiSR{x)a{x,s')] - fsoftjSRcris') (5.23) 

\Jxo / 0{a) 

for any s' in all expressions. The last two terms only cancel up to O [^\x=i Ax^ 
Equally, F^^^ should be read as 

F« a{s) = (f fisnix) a{x, s)) (5.24) 
for all higher-order terms. 



5.3 Implementation in WHIZARD 

The implementation of (Ttot.lSR+ in the Monte Carlo event generator WHIZARD is sim- 
ilar to the implementation of the fixed 0{a) correction described in Section 14.31 We 
use fisR (IF.14P as a user-defined structure function for each incoming beam with the 
scale Q = po cos 9 such that only collinear photons are described (cf. Section rF.2.3p . 
We then integrate over the xi, X2 dependent effective matrix element Meff (|5.ip ac- 
cording to Eq. ()5.2p . Note that, in contrast to the fixed order method, there is no 
(xi,X2) phase space slicing involved. Therefore, in the integration over the fisR soft 
region, the x-dependence of the matrix element is actually taken into account, so that in 
the code implementation additional contributions as given in (I5.23|) and ()5.24p appear. 
Note, however, that fsoft (|3-6p in M.eff still sets x = 1 throughout the soft region. The 
effects of this can be seen in the AEj dependence of (ytot,iSR{+) ^s discussed in Chapter El 

The 2^3 part of atotjSR+, 

I dx,fisR{xi;A9.,,^) [ dx2/isR(x2;A0^,^) X / dTs |Al2^3(s)P 
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Figure 5.3: Effective squared matrix element (arbitrary units) with the one-photon ISR 
part subtracted, for e~e~^ Xixt ^ function of the polar scattering 
angle 6 at ^/s = 1 TeV. Left figure: Helicity combination — h H — ; right 
figure: — h ++. Solid line: Born term; dashed: including virtual and soft 
contributions for AE^ = 10 GeV; dotted: same with AE-y = 0.5 GeV. The 
collinear cutoff is fixed at A6^ = 1°. 



is again generated in a separate run of WHIZARD , applying the soft and collinear cuts 
AEy, A9y and convoluted with the structure function fisR for each incoming beam. 
The scale is set to Q = po cos 6 as given above. As in the fixed order method, both 
runs can be combined for integration and event simulation. 



5.4 Difference to fixed order method 
5.4.1 Negative weights 

The resummation approach does eliminate the problem of negative weights: shifting 
the energy cutoff below the experimental resolution, such that photons are explicitly 
generated whenever they can be resolved, the subtracted effective squared 2 — > 2 ma- 
trix element is still positive semidefinite in the whole phase space. After resummation, 
the only potentially remaining source of negative event weights is the soft-noncollinear 
region. Negative weights are absent as long as 

0(l) X -In^^^WA^^) < 1, 

vr s 

where the 0(1) prefactor depends on the specific process. For the cutoff and parameter 
ranges we are considering here, this condition is fulfilled. 

Figure [5^31 shows \Ai(,ff\'^ for two different AE^ cuts. We see that the positivity condition 
is fulfilled. Since neither the inclusion of the ISR structure function nor the addition of 
the 2^3 part introduces further sources of negative weights, unweighting of generated 
events is now possible, so this method allows for realistic simulation at NLO. 
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5.4.2 Higher-order photon contributions 

As the fixed order method discussed in Chapter H] does not include higher-order contri- 
butions, the resummation method differs in all higher-order terms discussed in Section 
15.21 Below, we give the analytic expression for these differences. In addition, we discuss 
the difference between the resummation method and the inclusion of ISR as discussed 
in Section 13.71 Throughout the section, we neglect the contributions coming from the 
respective 2^3 integrations, i.e. processes where the last photon is hard and non- 
collinear. Here, differences between the methods are obvious from Section 15.21 



Radiation off one particle 

In this case, we again consider -A^g^j (|5.3|) and ct^^^jsr (|5.4|) . i.e. only subtract one 
factor of f so ft, ISR- We consider the soft and the hard region separately and again assume 
that, in the matrix element, x can be set to 1 in the soft region (j5.2ip . 

In the soft sector, we obtain 

A.oft = / dxflSK{x)\M^^^Wx,s) -(l + Aoft)|AlBorn|'(s)-2i?e(7WBornA^l_loop)(s) 

JXQ 

= (/soft.ISR (/soft - /soft,ISR) + F^"^^) \MBoTn?{s) + 2 Re (TWeorn Ali_loop)(s) /soft,ISR 

+ 0{a^). (5.25) 

In the hard sector, we obtain 

Ahard = [' dxfisR{x)\M'^Ji\Hx,s) - H dx T (x) \M^,j\x , s) 

xo 

dx {{{f^ix) {fsoft - fsoftJSR) + F2{x)) \MBorn\'^ix,s) 



fXO 

+ f+{x)2Re{MBo.nMl_,,,^{x,s)] -, / dxr{x)\M\ 

Jo 

+ 0{a^). (5.26) 
These results correspond exactly to the O(a^) contributions discussed in Section [Ql 

Radiation off two particles 

We integrate now over |A^e//|^ (15. ip according to (15. 2|) . 
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In the soft-soft sector, we obtain 
"1 ri 



^ss = dxi dX2 fisnixi) flSR{x2)\Meff\'^{xi,X2,s) 

J XO J XO 

-{l + fsoft\M Born 
= |2 {jsoft,ISR ifsoft - fsoftJSR) + F^"^^^ - fsoftJSR^ \M Born\'^ (s) 
+4/so/t,/5i?«e(7WBornA^lloop)(5) + 

= 2A,,ft-flf,^jsR\MBorn\\s) + 0{a^). (5.27) 



In the hard-soft section, we have 

/■I f^o ^ 

Ash = dxi dx2 fisnixi) flSRix2)\Meff\'^ixi,X2,s) 

Jxo JO 

rxo 



i-xo ^ 

- dX2Y. nX2)\M''Born\\x2.s) 

JO ^ = ± 

rxo 

= / dX2 {[f^ix2) ifsoft - '2 fsoftJSR) + F2{X2)] \M Bornl"^ ix2, s) 

Jo 

rxo 

+ /+(x2)2i?e(7WBornA^lloop)(a;2,s)} - / dX2f~{x2)\MBorn\'^{x2,s) + 0{a^) 

Jo 

= Ahard + 0{a^). (5.28) 
In the hard-hard sector, 

rxo rxo 



rxo rxo , 

Am = / dxi dX2 flSR{xi) flSR{x2)\Meff\'^{xi,X2,s) 

Jo Jo 

rxo rxo 

= / dxi dx2Fi{xi)Fi{x2)\MBornWxi,X2,s) + Oia^). (5.29) 
Jo Jo 



We see that, in general, nearly all second order contributions for the simultaneous photon 
radiation of e"*" and e~ can already be obtained from the second order terms for only one 
radiating particle. This only works in the approximation (j5.2ip . as this approximation 
does not distinguish the emitting particles for soft photons. We will make use of this for 
a rough higher-order estimation in Section I6.2.2[ 



Resumming only the Born matrix element 

We now consider the difference between (Jtot,iSR+ (|5.19|) and o"tat,/s'i?(f,) (|3.16|) : in the 
latter, only the Born cross section is convoluted with fisR- The leading and first order 
terms, which are already contained in the NLO corrected matrix element, are then sub- 
tracted to avoid double counting (cf. Eq. (I3.16P ). Up to 0(a), crtot,isR{b) corresponds 
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to the fixed order result atot (|3.1ip . As in the last section, we only consider O(a^) con- 
tributions. 

In the soft-soft region, the second order effects from convoluting the Born cross sec- 
tion are given by 

{flft,ISR + ^F^'^)\MBorn\Hs) 

and Ass (|5.27p becomes 

A'ss = 2 {fsoftJSR {fsoft - 2 fsoftJSR)) \MBorn\\s) + 4 fsoftJSR Re (MBorn M*_i„op)(s) 

+ 0{a^). 

In order to understand this expression, we again have to decompose fsoftjSR and fsoft 
into parts radiating form the same or two different particles and include a decomposition 
as (|5.1ip . obtaining 

2 UsoftJSR {fsoft — '2 fsoft,ISR)) = 

(^fW AFW) (xi) + (f(i) AFW) (x2) + F«(xi)AfW(2;2) + F«(x2) A 

This corresponds to the substitution of the fisR,soft by fsoft for the last soft-collinear 
photons in (Ttot,iSR+ and the substitution of the interference term description for virtual 
collinear last photons. In addition, 4 /so/t^/^R i?e (A^Born -A4*_[qqp)(s) gives all contribu- 
tions of (Trest, i-c. virtual contributions untouched by the subtraction, combined with 
the radiation of a soft-collinear photon from either incoming particle. 

In the hard-soft region, the 0{a^) contributions to crtot,iSR(b) are given by 

/ dX2 [f^{x2) fsoftJSR + F2{X2)] \MBorn\'^{x2,s) 

Jo 

and Ash (!5.28p becomes 

A's/j = / dx2 f'^{x2) {fsoft -'^fsoft,! SR)\MBorn\'^{x2,s) 

Jo 

f-XO 

+ /+(x2)2i?e (MBorn Mi*i„„p)(x2,s)} " / dX2 f {X2) \MBorn\'^ {X2, s) + 0{a^). 

Jo 

The hard-soft part of this is equivalent to 

F2{X2) ((AF«(xi) + AF«(X2)) . 

This corresponds again to the substitution of the fisR,soft by fsoft for the last soft- 
collinear photon (virtual last photon and interference term). 

In the hard-hard sector, the contribution to C(ot,/5i?,(6) up to 0{a^) is given by 
([529]). 
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In summary, we can say that the difference between the 2 — > 2 integration regions of 
(^tot,iSR{+) and atot,iSR{b) are the substitution of the fisR,soft by fsoft for soft-cohinear 
and the interference term description for virtual colhnear photons in (ytot,iSR{+)- In 
addition, terms of the form 



/ dxi / dX2flSR{xi) flSR{x2)c'rest{xi, X2, s) - drestis) (5.30) 
JO JO 

are included, where a rest is defined as in Eq. ()5.16p . These are the dominant contribu- 
tions, as all other differences between atot,iSR and crtot,iSR{b) are proportional to AF^^^ 
and therefore subleading. The O(a^) contribution to ()5.30p is given by 

2 dx f^{x)arest{x,s) + flSR,soft^rest{s)^ ■ (5.31) 

For e~^e~ — > xt Xi ^ these contributions are in the % regime; cf. Section [6. 2[ The sign 
of the correction term depends on the specific process (note that fsoft ISR < 0, while 
/+(x) > 0). 
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6 Results 



6.1 Cut dependencies 
6.1.1 Energy cutoff dependence 

In Fig. I6.1l we compare the numerical results obtained using the semianalytic fixed-order 
calculation with our Monte-Carlo integration in the fixed-order and in the resummation 
schemes, respectively. Throughout this section, we set the process energy to y/s = 1 TeV 
and refer to the SUSY parameter point SPSla'. All 2 — > 2 and 2 — > 3 contributions are 
included, so the results would be cutoff-independent if there were no approximations 
involved. For the discussion of the AE'^ dependence of atot, we refer to Section l4.6.1i 

The fully resummed result crtotjSR+ shows an increase of about 5 %o of the total cross 
section with respect to the fixed-order result which stays roughly constant until AE.y > 
10 GeV where the soft approximation breaks down. This increase is a real effect; it is 
due to higher-order photon radiation that is absent from the fixed-order calculation. 
We now discuss the reshuffling of contributions in the overlap region of the soft-collinear 
and hard-collinear (ISR) descriptions for second order contributions. If we raise AE^, 
photons that have been hard now become soft. In the case of photons radiated from two 
different external particles, a soft-collinear photon which has been described by fisR is 
now described by fsoft according to the soft approximation (|3.6p . For the case of two 
photons radiated off the same particle, we have to distinguish whether the two photons 
are A:^ -ordered or not. If they are, the description again changes from fisR to fsoft- If 
there is no -ordering, either the first or the second radiated photon can change from 
hard to soft (or both). For the first photon, it is a smooth transition where only the 
last two terms of Eq. (I5.20|) are involved. If the second photon changes to soft, new 
contributions of the form AF^^^ (I5.17P appear. 

The effects of reshuffling the photon descriptions are, up to AE^ = 10 GeV, smaller 
than 2%o. If we only consider the dominant k± contributions, it would be preferable to 
raise AE^ as long as the 0{a) effects from the soft approximation are negligible. In our 
case, these effects reach the 2%o level for AE-y = 1 GeV (cf. Section [1. 6. ip . 
Finally, we compare the difference between atotjSR and (7tot,iSR+- For the latter, the 
emission of a last hard non-collinear photon is also combined with multiple photon 
radiation described by fisR- We observe that for AE^ > 1 GeV these higher-order 
contributions are caught by ISR resummation of the 2^2 part. For a lower cutoff, 
they have to be explicitly included in the 2 — > 3 part. 
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Figure 6.1: Total cross section dependence on the energy cutoff IS.E^ using different cal- 
culational metfiods: 'sa' (red, dotted) = fixed-order semianalytic result using 
FeynArts/FormCalc; 'fix' (green, dashed) = fixed-order Monte-Carlo result 
(^tot (|3.1ip using WHIZARD; 'res' (blue, long-dashed) = ISR-resummed Monte- 
Carlo result crtot,lSR+ (|5.19p using WHIZARD; (magenta, dash-dots) = same 
but resummation applied only to the 2 — > 2 part (j5.18p . Statistical Monte- 
Carlo integration errors are shown. For the Monte-Carlo results, the collinear 
cutoff has been fixed to A9.y = 1°. The Born cross section is indicated by 
the dotted horizontal line. 
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Figure 6.2: Relative dependence of the cross section on the coUinear cutoff using 
different calculational methods: sa (red, dotted) = fixed-order semianalytic 
result using FeynArts/FormCalcwith respective integration errors; fix (green, 
dashed) = fixed-order Monte-Carlo result otot (|3.1ip using WHIZARD; res (blue, 
dash-dots) = third-order completely resummed Monte-Carlo result (Jtot,iSR+ 
(I5.19P using WHIZARD; (magenta, long dash-dotted) = third-order resummed 
Monte-Carlo result (JtotjSR (|5.18|) using WHIZARD; Statistical Monte-Carlo 
integration errors are shown. The soft cutoff has been fixed to /S.E^ = 
10 GeV. All results are scaled to a Born = 1- 



6.1.2 Collinear cutoff dependence 

In 0{a), the collinear angle dependence tests the validity regime of the collinear ap- 
proximation used in Eq. (j3.9p . This was already discussed in Section 14.6.11 where we 
found the approximation is valid for < 1°. The difference between the fixed order 
result atot and the fully resummed result aisR,+ are again the higher-order contributions 
discussed in Section 15. 2[ 

Figure [62] shows the higher-order effects associated with raising the collinear cutoff A0^. 
For = 0.1° (1°), the difference between atot and (Jtot,iSR+ is 4%o(7%o). Raising 

A0^ shuffles photons from a non-collinear to a collinear description. First, this opens up 
phase space for all photons described by fisR- For the radiation off the same particle, 
raising the collinear cutoff also improves the description of ordered photons if the 
sum of their transverse momenta is now lower than the new transverse cutoff associated 
with A0^. They are contained in the second order ISR contributions F^'^\ F2 and switch 
from a description oc A F^-^^ to terms of the form (j5.20p . where the second soft (virtual) 
collinear photon is given by fsoft (the interference term). 

Unfortunately, raising the collinear cutoff also enhances the phase space of contributions 
oc AF(^) which can appear when multiple (collinear) photons which are not k± ordered 
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Figure 6.3: (JBom-, o'tot, criSR,+ as a function of ^/s. All results in fb 

are emitted from the same particle; cf. Section 15.21 Similarly, raising the cutoff can 
worsen the description for photons that lie in the soft regime near the limit of the soft- 
collinear regime and change into the latter. For non-fe^-ordered photons, the description 
switches from fsoft (interference term) to contributions oc AF^^\ However, these effects 
are subleading. 

For the leading-log contributions, raising the cutoffs gives a better description as long as 
the collinear approximation is valid. From Figure [6^ we see that this is a 3 — 4%o effect. 
The difference between atotjSR and crtotjSR+ is similarly small (about 1 %o) in the va- 
lidity regime of the soft approximation. The size of these contributions therefore mainly 
depends on the choice of the soft cut AE^ (cf. Section [6.1.ip . For the AE^ = lOGeV 
as above, all leading contributions are already summed up in the ISR part of the 2 — > 2 
integration. 

6.2 Total cross section 

6.2.1 Leading and first order results 

In leading and next-to-leading order, the fixed order result atot (|3.1ip and the resum- 
mation result (ytot,iSR(+) (|5.18p . (15.190 agree. All differences between the results are 
therefore due to higher-order effects. For the discussion of first order behavior of the 
total cross section, we therefore refer to the Sections 13.61 and 14. 6[ 

In Figures 16.31 and 16.41 we display the LO result together with the NLO results for the 
fixed-order and resummed approach. Near the cross-section maximum, the relative cor- 
rection in the resummed approach is about —5.5%, approaching —2% at y/s = 1 TeV. 
Near threshold and at asymptotic energies, the relative NLO correction gets up to —20 %. 

A short remark concerns the neglected helicity-flipping terms f~{x) (jF.Sh in the hard 
part of fisR which are in principle ©(a), but do not contain collinear logarithms; cf. 
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Figure 6.4: relative difference of atot (and crfSR+) to (JBom, resulting from NLO (NLO 
and higher-order) corrections 

Section IF. 31 Figure 16.51 shows the relative magnitude of these terms for unpolarized 
initial particles with respect to the Born cross section. The contribution is in the per 
mille regime. 

6.2.2 Higher-order effects 

Figure 16.61 shows the differences between atotjSR-, o'tot,iSR+j (^tot,iSR{h) aiid the fixed 
order result atot p.lip . atotjSR, CFtot,iSR+, and atot,iSR(b) all contain different kinds of 
higher-order contributions (cf. Section [5.4.20 . which we will now consider in more detail. 
The largest effects clearly result from the higher-order contributions in the transition 
from the fixed order result atot (13. Hp to the resummation result atotjSR (I5.18P discussed 
in Section 15.4.21 i.e. the inclusion of multiple photon radiation off the Born as well as 
the interference term in the 2 — > 2 integration region. They are up to a few percent for 
small center of mass energies. Including the convoluting of the 2 — > 3 integration region 
of atotjSR with fisRj i-e. including also diagrams where the last radiated photon is hard 
and non-collinear as in atot,iSR+ (!5.19p . significantly reduces these effects. These extra 
terms cannot be neglected. However, it was shown in Section [6. II that the magnitude of 
these corrections strongly depends on the soft cut Ai?^. If the energy cutoff is high 
enough, the dominant contributions are included in the convolution of the effective 
matrix element in the 2^2 region with fisR] cf. Figure 16.11 Apparently, the cut 
value IS.E^/y/s = 5 • 10~^ chosen for Figure is not sufficiently high for low ^/s 
values to include all dominant contributions. 

Finally, we see that the difference between atot,iSR{h) and atot,iSR+-, which is mainly given 
by the terms oc / arest (I5.30P . is similarly in the % range. Shuffling photon descriptions, 
in contrast, is a %c effect; cf. Section [6.11 The change of ISR order in the resummation 
method only concerns the hard regime and clearly shows no significant effect. 
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Figure 6.5: Contribution of the part of atot resulting from the hehcity-flipping part /~ 
()F.5p of the hard-cohinear approximation with respect to (JBom- 1 c errors 
from WHIZARD integration are shown 
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Figure 6.6: Relative higher-order effects for different methods: (magenta, long dash dot- 
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Third (and higher) order estimation 

Eqs. (|5.27p . (j5.28p . and (j5.29p . show that we can in principle do a 0{a'^) and higher- 
order estimation when using the results from radiation off one particle (|5.25p . 
()5.26p to predict the O(a^) effects for the radiation off two particles. Analytically, we 
have 

A,, + 2Ash = 2 {As + Ah) - fsoft,ISR\MBorn\^ + 0{a^). (6.1) 

Deviations from the "estimate" on the right hand side can give us an idea of the order 
of magnitude of the higher-order corrections. Note, however, that numerically the cal- 
culation of these terms are very involved and limited by computer precision; similarly, 
errors resulting from setting x = 1 in the matrix element the soft regime are neglected. 
This estimate also does not predict the actual sign of the higher-order corrections. It 
tests the numerical difference between expressions which are analytically equivalent up 
to 0(q2). 

Figure [UTTl shows the estimated and actual 0{a'^) contributions, here for the convolution 
of \M-Born\'^ as in /5/j(b) ()3.16p only where similar considerations can be done. We 
see that the deviation from the prediction is already in the 0.5 — 1% regime. Note that 
for a comparison we also added points where 

f'XQ f'XQ 

A^r = dxi dX2flSR{xi)flSR{x2)\MBorn\^ 

Jo Jo 

was replaced by Ahh (15.290 . These two expressions differ in third and higher-order terms 
involving the radiation of two and more hard-collinear photons off the virtual (+ soft- 
collinear) NLO term. They are of similar magnitude. 

FigureESlshows the magnitude of the 0{a^) errors in Eq. (16. ip relative to the Born cross 
section and the equivalent contribution from cr^^^ jg^K^f,-^ in Figure [H771 We see that, for 
^/s < 600 GeV, the third order effects are of similar order for both methods. For larger 
cm energies, however, the estimate from atotjSR deviates from the Born-related estimate 
and goes up to nearly 2%. This is related to second and higher-order effects resulting 
from the convoluting of the virtual part with fisR for both beams. The importance of 
these contributions was already shown in O(a^) in Figure [6^61 as well as the difference 
between A^^™ and A^h in Figure 16.71 Summarizing, we can say that according to 
the estimation done in Eq. (16. ip . even third and higher-orders are important if the ILC 
precision of a few %o is reached. 

6.3 Event Generation 

As for the total cross section, events generated using M.eff (14.40 and Meff (15. ip are 
equivalent up to and including 0{a). For the discussion of the NLO effects, we therefore 
refer to Section 14.61 Differences between the two methods originate from higher-order 
contributions. 

Figure 16.91 shows the binned distribution of the chargino production angle as in Figures 
14.81 and 14.91 together with the events generated using the resummation method. For a 
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Figure 6.7: Second and higher-order effects; estimate from 1 photon (blue, dashed) com- 
pared to actual result (red, solid) . Born only 
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Figure 6.9: Polar scattering angle distribution for an integrated luminosity of 1 ab^ at 
^/s = 1 TeV. Left: total number of events per bin; right: difference w.r.t. 
the Born distribution. LO (black, dotted) = Born cross section without ISR; 
fix (green, dashed) = fixed-order approach; res (blue, full) = resummation 
approach. Cutoffs: AE^ = 3 GeV and A6^ = 1°. 



kinematically more accurate description of the NLO contributions, decreasing the cutoffs 
would be preferred, but choosing lower values would invalidate the fixed-order approach 
for the comparison. 

Figures I6.1UI and 16.111 show the difference in the number of events between the angular 
distribution for fixed order and the resummation method as well as a 1 cr effect from the 
Born cross section. The higher-order effects are not as striking. They are visible mostly 
in the central-to-forward region. For cos 9 > 0.5, they are statistically significant. Going 
to higher luminosities would of course improve the significance for even lower values of 
cosO. From Figure [6. Ill we can see that the integration over cos^ reproduces the ~ 6%o 
effect shown in Figure [6T] (compare also to Figure \6M with AE^ = 5 GeV; from Figure 
16. H we know that the difference between the resummation and the fixed order result 
stays roughly the same up to AE'^ = 10 GeV). 

Finally, we can test the effect of lowering A E^ such that we get into the critical soft cut 
region discussed in Section ITSl Figure [6.121 shows the difference between the two methods 
of adding photons for AE^ = 0.5 GeV, AO = 0.5°. We clearly see the effects of setting 
the negative effective squared matrix elements to zero where it becomes negative. The 
behaviour of |A1e//P for the (sub)dominant helicity combination was shown in Figures 
14.21 and 14.31 For cos^ > 0, the differences between the resummation method and the 
fixed-order method have a behaviour as in the higher energy cut case, cf. Figure 16.101 
They are smaller in magnitude. As in this region of phase space, |A^e//|^ ^ for the 
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Figure 6.10: Polar scattering angle dependence of difference between events resulting 
from completely resummed and fixed order method: N^es ~ ^ex- Standard 
deviation from Born (solid, red) is shown 




Figure 6.11: Same as Figure [6.101 normalized to Born results 
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Figure 6.12: Difference between complete resummation and fixed order metliod for ^/s = 
ITeV with AE^ = 0.5 GeV, AB^ = 0.5°. With the energy cut in the 
critical region (cf. Section effects of setting |7We//P to zero become 
apparent for cos 9 < 0. 



dominant helicity combination, this difference can be attributed to the lower angular 
cuts which decrease the difference between atot and crtot,iSR+ (cf. Section [6.1.2p . 
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7 Summary and Outlook 



As an e^e~ collider with a cm energy of 500 GeV (ITeV), the ILC will provide a clean 
environment for precision tests of physics in and beyond the Standard Model. In the 
chargino and neutralino sector of the MSSM, a small number of measurements already 
determines the supersymmetric parameters at the electroweak scale. The parameters 
at the SUSY breaking scale can be derived from renormalization group equations and 
eventually give information about the SUSY breaking mechanism. The clean experimen- 
tal environment of the ILC leads to small experimental errors which calls for matching 
accuracy of theoretical NLO and higher-order predictions. These effects need to be in- 
cluded in Monte Carlo event generators, the simulation tools that are actually used in 
the experimental analyses. 

Inclusion of NLO corrections in WHIZARD 

We have presented results obtained from implementing the NLO corrections for chargino 
pair production at the ILC into the Monte-Carlo event generator WHIZARD . On top of 
the genuine SUSY/electroweak corrections, we have considered several approaches of in- 
cluding photon radiation, where a strict fixed-order approach allows for comparison and 
consistency checks with published semianalytic results in the literature. However, this 
approach suffers from negative event weights in certain points of phase space. Extending 
approaches used in Monte Carlo generators for LEP analyses, we developed a method 
for including virtual-, soft- and hard-collinear resummation of photon radiation. It not 
just improves the numerical result, but actually is more straightforward to implement 
and does not suffer from negative event weights in or near the experimentally accessi- 
ble part of phase space. The generator accounts for all yet known higher-order effects, 
allows for small cutoffs, and explicitly generates photons where they can be resolved ex- 
perimentally. In addition, the generated event samples allow for quick analyses of NLO 
effects in angular distributions, correlations, and other quantities. We also interfaced 
WHIZARD with a (modified) FormCalc code which in principle allows for the inclusion 
of any other (NLO) process. 

Magnitude of next-to-leading and higher-order corrections 

For the mSugra point SPSla', the NLO corrections for chargino production and decay 

cross sections at the ILC arc in the percent range. Convolution of the virtual non- 
collinear part with the ISR structure function leads to additional percent effects for 
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large ^/s. These contributions are not included in the standard treatment of ISR resum- 
mation, where multiple photon emissions are combined with the Born process only. 

Minor next-to-leading and higher-order uncertainties for the cross sections stem from 
the use of approximations in the soft and collinear regime, reshuffling of photon descrip- 
tions, enhancement of phase space for higher-order photons, and neglecting of the non- 
logarithmic 0{a) photon radiation from the helicity flipping term in the hard-coUinear 
approximation. A rough estimate shows that even higher-order contributions described 
by the resummcd structure functions can become C'(%) and are therefore important in 
a thorough analysis of the process and dominant with respect to the non-leading helicity 
flip contribution. The latter is O{%o) for unpolarized cross sections; it can in principle 
be easily implemented in the resummation method. 

Cutoff dependencies 

In general, a careful analysis of the dependence on the technical cutoffs on photon energy 
and angle is necessary, as it reveals uncertainties related to higher-order radiation and 
breakdown of the soft or collinear approximation. The soft approximation leads to a 
deviation from the exact solution of ~ 5 %o for a ratio of the soft photon cut to the 
center of mass energy I^E^/y/s = 10~^. The collinear approximation breaks down for 
the collinear cutoff A ^ = 1°. We compared the corrections resulting from the use of the 
approximations to the magnitude of the NLO effects as well as higher-order corrections. 
For the precision reached at the ILC, both the soft and collinear approximations imply 
the use of low cuts. On the other hand, raising the cutoffs in principle enhances phase 
space for multiple photon emissions and reshuffles photon description from the leading log 
initial state radiation description to the soft approximation for soft and the exact virtual 
description for virtual photons at leading order. These reshuffling effects are in the %o 
regime, but smaller than the errors arising from raising the cutoffs for the first order 
approximations. This implies the use of low cutoffs for a more exact NLO description. 
The 2^2 process describing the 0{a) emission of hard non-collinear photons also 
needs to be convoluted with the ISR structure function, as these contributions also are 
in the %c regime. 

Cross sections and event generation 

With both the fixed order and the resummation method implemented in WHIZARD , we 
reproduce the semianalytic literature results for chargino production at NLO. The angu- 
lar event distributions reflect the size of the corrections for total cross sections discussed 
above. For a center of mass energy of ITeV and an integrated luminosity of lab~^, 
simulation results derived from both the NLO fixed order and the resummation method 
differ significantly from the Born distribution and each other. The corrections cannot be 
described by a constant proportionality factor K between the Born and the NLO result. 
The fixed and resummed distributions differ in higher-orders. This is most obvious in 



73 



the forward-scattering region, where the distributions differ by more than three (Born) 
standard deviations. 



Outlook 

The constructed generator should be regarded as a step towards a complete NLO sim- 
ulation of SUSY processes at the ILC. If charginos happen to be metastable, it already 
provides all necessary ingredients. Beam effects (beamstrahlung. beam energy spread, 
polarization) are available for simulation and can easily be included. However, charginos 
are metastable only for peculiar SUSY parameter points. If the theoretical prediction 
should match the experimental accuracy, a full description of the end products seen in 
experiment needs to include the chargino decay at NLO as well as non-factorizing contri- 
butions to the 2 ^ n process, e.g. in the double-pole approximation. These we have to 
match with off-shell and background effects, already available for simulation in WHIZARD. 
Furthermore, in the threshold region the Coulomb singularity calls for resummation, not 
yet accounted for in the program. These lines of improvement will be pursued in future 
work. 
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A Conventions, SUSY overview, MSSM 
Lagrangian 



There are various introductions to the formal construction of a = 1 SUSY Lagrangian, 
cf. [301 EH [2ni I139t [331 ED- We just cite the main results and refer to the literature for 
further details. 



A.l Conventions 

In the following, we define g'^jy by 

5m. = 9^" = diag(l, -1,-1,-1). 
The Dirac matrices obey 

{iti^lv} = '^9lMU l4x4- 

In the Weyl representation, they read 

■'^ - ( ) . <-) 

with 

= {l2,a') , = {l2,-a') , (A.2) 
and the being the Pauli matrices. In addition, we define 

4 4 
and 

The two-component Weyl spinors transforming under the (^,0) and (0, ^) representation 
of the Lorentz-group are given by 

^« (^,o),r (0,^). 

The index-raising and lowering matrices are given by 

^ ~ ~ V 1 

ea, = = ( J ) • (A.4) 



(A.3) 
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A four-component spinor is then given by 

The index structure of the a matrices (IA.2j) . (1A.3P is 
in accordance with Eqs. (jA.ip and (|A.5p . 

Multiphcation of Weyl spinors (as well as any Grassmann numbers) is defined by 



Together with Eq. ()A.4h . this leads to 



A. 2 Poincare and SUSY algebra 



We define the generators of the translation group and of the Lorentz group such 
that the infinitesimal transformations are 

x'fj, = Xfj, + ttfj, = x^-iax{P^)^ (translation), 

x'P = xP + ujP„ x" = xP uJ^,y{M'"')P^x'' (Lorentz transformation). 

The generators P^, M^jy then obey the Poincare algebra 

pX pfJ. 



0, 



Mi'^P^ = i{g''^P>' - gf^^P"), 
[Mi"',MP''] = iig^PMP" + gf^^M^P - gPPM"^ - g^^MPP). 

A supersymmetry transformation is an extension of the Poincare Algebra relating bosonic 
and fermionic degrees of freedom. We define Qa, Q" to be the supersymmetry generators 
in the two-dimensional (^,0), (0, ^) representations of the Lorentz group. The properties 
under transformations and closure of the algebra require (anti-) commutation rules which 
in four-component notation read 

[M'^^Q] = 0, {Q,Q} = 0, [M^^Q] = -^^p" Q, 

where 



76 



In two-component notation, this reads 

{Qa,Q/3} = = 0. (A.6) 

A. 3 Expansion in component fields 

A SUSY transformation for a field operator (j){x) is given by 

exp [i iCQa + Q")) (Pix) exp (-i iCQa + 

where ^ are Grassmann-valued variables. Infinitesimally, this becomes 

<Pix) ^ 0(x) + 6^4>ix) , S^^ix) = [i{CQ + ^Q),Hx)] . 

Starting with a scalar field (j){x) and requiring the closure of the algebra, we obtain the 
transformation rules 

5^F = iV2ia^'d^^, (A.7) 
where "0 and F are a spinor/ auxiliary field with the dimensions |/ 2 respectively. 

From (|X6]) . we can determine the action of a SUSY transformation generator 

G(a^ e, = exp (i (e Q + - a^^^^.)) 

on a superfield S{x^,9,9): 

S{x^,e,e) ^ exp(i(eQ + eO-«^^'/.)) 5(2;'^, e,^). 

S{x^,6,6) can be expanded in components proportional to 9, 9, 99,...; the generators 
are then given by 



d_ 
d_ 



Fermionic derivatives which anticommute with the Q, Q are then given by 

d . ^ _d_ 

' dt 



Do, = —+ia>'9d. , = - i9a^'d., 
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which obey 

All other anticommutators vanish. 



Defining a chiral superfield <I> by requiring 

^ = 0, 

we obtain 

e,e) = + ^Oii^ + eeF + idf, (poa^'e - eod^ii: a^^e - \df, a'^ < 

v 2 4 



with the component fields and their transformation rules given by Eq. ()A.7p . 
The antichiral field is given by 



V 2 4 



which obeys 



A vector superfield V is defined by requiring 

V = 

This leads to 

V{x,e,e) = C{x)+i9x-i9x + ^iG&[M{x)+N{x)]-^i99[M{x) -iN{x)] 



9a''9Vf,{x) +i999 



i999 



+ - 



1 



D - -a,9^c 



with C, M, A^, D being real scalar fields, Xi ^ Weyl spinor fields, and a real vector 
field. Note that there is a gauge invariance with respect to 



V V + 



which helps to reduce V to 



Vwz{x,9,9) = 9a^9V^{x)+i999\{x)-i999\{x) + 



D 
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($i$2), ($i$2^'3), {D^D^V) chiral fields 
($|^>i), i [$1 - vector fields 

Table A.l: superfield character of different superfield combinations 

in the Wess Zumino (WZ) gauge. The fields V^, A, D transform under SUSY transfor- 
mation as 

with the field strength V^^, = d^Vi, — d^V^. We see that supersymmetry is explicitly 
broken in the WZ gauge; however, it can be restored by a proper gauge transformation. 



Combinations of (anti)chiral and vector superfields which are again (anti-)chiral and 
vector superfields are given in table lATl A more extensive list can be found in [32j. 

A. 4 General gauge-invariant Lagrangian and superpotential 

For the construction of a Lagrangian, we have to find components of superfields which 
are super symmetric. It is easy to see that these are given by the 99 {99) components 
of the $($^^) (F-terms) and the 9999 components of V (D-term) respectively, as all of 
these transform up to a total derivative. Lagrangians are therefore constructed using 
the F-terms of chiral and the D-terms of vector superfields and the field combinations 
as given in table lA.ll 



Gauge invariance 

For a general U{1) invariant Lagrangian, we define the field-strength superfield 

9) = AiX^{y)-[A5jD{y) + 2i{a^^anJV^,{y)]9p + ^99ard^X^, 

where = + i9a^9. This contains a spinor with the dimension | and the field 
strength V^,y with dimension 2 as well as the auxiliary field D. The pure gauge term is 
then given by the F-term of (WaW^ + he): 

Cv =^iW''W^ + WjV'^)F = -^V^'''V^,+ iXa^'d^X + ^D^ (A.8) 

The interaction between matter and gauge bosons is then given by the D-term of the 
vector superfield ^: 

C^ = [<^^e^^]D = {D^(l))\D^'(l))+iiPa>'Dl^ + F^F + iV2g{(j)'^Xilj-^X(l)) + g(l)'^Dcl), 

(A.9) 
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where = + i gV^. Then, the total gauge-invariant Lagrangian is 

For the construction of QED with massive fermions, we have to include left- and righthanded 
component fields; this can be done by defining the combinations 

1 1 

$+ = ^ (<I>1 1^2) , = ^ ($1 - i^2)- 

Their transformations are then given by 

$+ ^ exp(-2iA)^>+ , ^ exp(2iA)$_ ,V ^ V + i(A-A"^). 
In addition to , £$_ , this allows the inclusion of mass generating terms of the form 



so the total gauge-invariant Lagrangian including the coupling of massive particles to 
the gauge bosons is 

C^inv — C^V "I" + ^<^- + -^mass- 

This procedure can be extended to any (non)abelian gauge group. We have to add terms 
in the form of Eqs. (jA.Sp . (|A.9P for every new gauge group and modify accordingly. 
In the non-abelian case, the vector superfield is given by 

where f""^"^ are the structure constants of the non-abelian gauge group. Its transformation 
is given by 

Superpotential 

The superpotential defines the possible self- interactions terms of chiral fields. Renor- 
malizability requirements determine the most general SUSY-preserving potential 

i i,j i,j,k 

Elimination of auxiliary fields 

The unphysical auxiliary fields Fi,D are generally eliminated by requiring 

dC _ d£ _ 

dD ~ m ~ ' 

which leads to 

= -§T.^D = -Y,49<P^- (A.IO) 
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gauge group 


non-singlet SM matter fields 


U{1) (Y) 

SU{2) 

SU{2,) 


UL,dL {\), (-^), UR (-§), dn (i), Ir (1), H (i) 

UL,c, dL,c, Ur^c, dR^c 



Table A. 2: Standard Model matter content and gauge-group transformation properties. 

u = {u,c,t), d = {d,s,b), I = {e,iJ,,T) run over the respective families and c 
is the color index. Y is the hypercharge. Transformation properties of gauge 
fields are not listed. 

A. 5 MSSM field content and superpotential 

As the MSSM is the minimal supersymmetric extension of the Standard Model, it is 
designed to reproduce the SM particle content including as few extra fields as possible. 
We therefore need 

• Vector fields as in (1X8]) for the SU{2) x C/(l) electroweak as weh as the SU{3) 
strong gauge groups 

• gauge-particle interactions as in (|A.9h for the Standard Model matter content as 
hsted in table IA.21 

• a superpotential Jmssm giving masses to the (s)fermions. This requires a modified 
Higgs sector: instead of the SM Higgs, the two Higgs doublets given by Eq. (IA.12P 
are introduced. They carry the hypercharges Y = {Hu) and Y = — ^ {Hd}- 

The minimal superpotential reproducing the Yukawa-terms of the Standard Model and 
preserving supersymmetry is 

fMSSM{<P) = uyuQHu - dydQHd - eyeLHd + fiHuHd, (A.ll) 

with 

and u, d, e denote the righthanded SU(2) quark and lepton singlets and Q, L the left- 
handed SU(2) quark doublets, yi are the respective Yukawa matrices. 

A. 6 Symmetry breaking and Lagrangian of the MSSM 

SUSY and gauge breaking 

Supersymmetry breaking terms need to preserve the cancellation of the quadratic di- 
vergences of the corrections to the Higgs mass. The most general terms fulfilling this 
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requirement are 

SUSY -breaking = + mlluij^ + mjIdRp + mf\!L\'^+mj\eR\'^ 

+ (^XeAeHJicr + \d AoHdqiUR + \u Au HuqhdR + BfiHuHd + h.c^ 
+ mljHul^ + mljHdl^ + \ (^M^B B + M2WW + M^~g ~g + h.c^ . 

(A.13) 

Here, m?, m|, m^, m?, m| are general hermitian 3x3 matrices and \e Ae, Xd Ad, Xjj Ajj 

general 3x3 matrices, \qL\^ = {u^j^ ul + d\ di) (same for II) and the multiplication of 
two SU (2) doublets is given by 

Di D2 = e'Wi^i D2,j = Z)i,i D2,2 - ^2,1 Di^2- (A. 14) 

As in the SM, the electroweak gauge symmetry is broken because the Higgs doublet 
aquire non-zero VEVs. However, the specific form of the MSSM Higgs potential requires 
both II in /mssm and B 12 in CsusY -breaking to be non-zero. 

MSSM Lagrangian 

The complete MSSM Lagrangian is then given by 

C^MSSM = -\ V^rVg.,,, + t \g. d^\, + - \ Vl^^^Vl^, + d,Xl + ^{Dlf 

+ {D^cI)i)\D^'^i) + i a'^Dli^j + f\Fi + i V2gk {^T'^X'^t/Ji - ^iX^T^<t)i) + gk^DKl^i 

I 2|~ |2 I 2\~ |2 I 2|j |2 I 217' |2 , 2|~ |2 

+ mq\qL\ + mf^\uR\ +m^\dR\ + m^\lL\ + m^\eR\ 

+ (^Xe Ae HjieR + Xd AoHdqiUR + Xu Au Huqid-R + BfiHuHd + h.c^ 

+ mj,jHu\^ + ml^ \Hd\^ + ^ (^Ah B B + M2WW + Ms g g + h.c}j 
+ (uyuQHu - dydQHd - eyeLHd + nHuHd)gQ , 

where i = 1, 2, 3 is the SU (2) and a = 1, 2, 8 the SU (3) group index; the index / in 
the gauge-matter interaction part ()A.9p runs over all superfields corresponding to the SM 
particles listed in Table IAT2] and the Higgs fields (|A.12p . where the index k symbolizes 
that all respective gauge-groups as given in Table IA.2I have to be taken into account. 
The field content from the superpotential is given by its (66) component and can be 
derived using Eq. (TOll) and 

^i^j\{ee) = -i>i^j, 
^i^j^k\(8e) = -{ipiipjcpk + cyclic terms). 
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bosonic 



fermionic, spin = ^ 



eL,R, fiL,R TL,R (sleptons, spin = 0) 

UL,R, dL,R CL,R sl,r, iL,R bL,R (squarks, spin = 0) 

B, W^, g°- (gauge bosons, spin 

H^, H^, H;^, (Higgs bosons, spin 



UL,R,dL,RCL,RSL,RtL,R. bL,R 

1) B, W\ g" 

0) H+,HtH^,Hl> 



eL,R, fJ'L,RTL,R 



(leptons) 
(quarks) 
(gauginos) 
(Higgsinos) 



Table A.3: field content in the MSSM 



Finally, D and F terms are eliminated using Eq. (jA.lOp . 

The field content of the MSSM is then given in Table Physically observed particles 
correspond to superpositions of these fields carrying the same quantum numbers which 
are rotated into mass-eigenstates. In addition, the two-component Weyl-spinors A, V 
have to be combined to obtain four-component Dirac- or Majorana spinors. For the 
gaugino/ Higgsino sector, this is done in Appendix [Bj For a more complete treatment, 
we refer to e.g. ^31)1 [33l l3i] . 
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B Chargino and Neutralino Sector of the 
MSSM 



Here, we sketch the derivation of the Feynman rules for any vertices including charginos 
and neutralinos. We hereby closely follow We point to differences in conventions if 
necessary. 



B.l Chargino mass eigenstates 

In the MSSM, the EW symmetry is spontaneously broken by non-zero VEVs of the two 
neutral components of the Higgs doublets: 

<Hl>>= VI, <H^>= V2. 

tan (3 is defined as the ratio of these VEv£|: 

tan p = — 

Vl 

Taking this into account in the chargino-chargino-Higgs couplings and using 
we obtain the following mass terms in the Lagrangian: 



{M2W+W~ - tiH+H^) + h.c. - igH^W-vi + igv\H;;^W- - igH^W+V2 + igv2H+W' 
We now define and ^'^ by 



and obtain 



as mass terms of the Lagrangian with 

X = 



M2 gv2 

gvi n 



^Note that the definition of tanS„ in [30] is given by tan0„ = (tan/3) ^. 
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Using the relation 



we can rewrite this as 



2 1 2/ 2 I 2\ 



X=i \/2m^sin/3 > 

* \/2mu, cos /3 /i ' ' 



This matrix can now be diagonahzed using 
with 

cos ^± sin ^± 



' - sm ^-i- cos *± 
The angles cos^'± and sin<^-|- obey [76l 178] : 

cos = - [Af| - ± cos 2/3] /Ac, 



sin2^>± = -2mi4/ VM| + |//|2 =f(M|-|/x|2) cos 2/3 + 2 M2I/XI sin 2/3/Ac, 



Ac = -^(M| - |/i|2)2 + 4m2^ cos2 2/3 + 4m2^(M| + |^|2) + 8m2,M2|Ai| sin 2/3. 
The two-component charged mass-eigenstates a-re then given by 

and the four-component Dirac spinors xt by 

,7 = (I) or XT = (I). 

depending on the particle/ antiparticle choice. 

Note that there are different notations for these matrices in the literature which have to 
be taken into account when comparing results. We here give a short overview: 
We start with the form 

X j U- 



with as by Eq. dHUand X as by Eq. (IM . We now define the matrices U, V such 
that 

X- = U^- ; x+ = 
and see that X is diagonahzed using 

u* xv\ 

Depending on the parameterization of U and V , m~± can be larger or smaller than m~±. 
Table ET] gives an overview of the different conventions used in the literature. 
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HK[30j 


V 


u 


u* XV^ 


here 


^ + 


Ul 




LH|129| 


V 


u 


u* xv^ 


Choi[78] 


Ur 


Ul 


uixul 


FC[80] 


V 


u 


u*xv^ 



Table B.l: Different conventions for chargino-diagonalization matrices in the literature 



B.2 Neutralino Mass Eigenstates 

From Cmssm, the neutralino mass terms are given by 



+^v^H^,W' - '-j=v,H',B ] + h.c. 



We define by 

\ m J 

with the four-component Majorana-spinors 



wl 



mi 



We can then rewrite the mass part of the Lagrangian as 



-iB\ ~ 



w 



--^l^ M-^l + h.c, 



with 



M 



/ 


Ml 





g' 


g' 


\ 







M2 





















V 











/ 



Using 



1/ 2 I /2\ 2 u 

7,\9 +9 )v ; — = - 



g_ , sp_ 

9 ' cp 



V2 2 

— ; V 



Hi 
Hi 



2 I 2 
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we obtain 



M 



/ Ml 




'w 







-w 



-mzcpSu, mzSfsSw \ 

mzcpcw -mzspcy, 

-II 

-^^ / 



This can be diagonalized using a unitary matrix such that 



diag 



N* MNl 




and the four-component spinors by 




When taking N purely real, neutralino masses can come out negative. In this case, the 
physical field is given by 75 x''- Note that this has to be taken into account when deriving 
the Feynman rules for the physical particles from C. 

An analytic expression for the diagonalization matrix can be found in the literature, 
cf. e.g. [HO] and [TIT] . 

B.3 Chargino-Gauge-boson couplings 

We will here and in the following switch between two-component Weyl spinors and four- 
component Dirac- and Majorana-spinors; a good introduction can e.g. be found in [139J. 
The part of Cmssm describing chargino-chargino-gauge boson coupling is given by 



-gH+am'^Ht 



--g'Hta^'B^H+ - gW+a^^Wj^W, 



u 



■+ 



+ gW-a^W'^W- 



Let us define the as the particle. If we use 




the Lagrangian is given by 




(B.3) 
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We can now split both W and H into left- and right-handed parts as 

and same with and H. We then obtain the following respective relations between 
the left- and righthanded parts of xi,2 and H, W^: 



X2L 



\Hl) ' \X2r) \Hr 



The whole coupling term of the Lagrangian containing left-handed parts now reads 



( g-i^'W^^ \ (Wl 

\hJ V ^(g^>^wji + g'rB^) J\HlJ- 

The same relation (with L <^ R) holds for the right-handed part. 

If we define the negatively charged charginos to be the particles, the derivation is similar; 
however, we then have to define 

\{-iW+)) \Hi 

The terms in the Lagrangian describing the gauge-boson couplings are the same up to a 
sign change. The transformation rules for the right- and left-handed parts are given by 

XI l\ ^jj-y(^L\ . fxiR\ ^ ^^^.JWr 



kX2LJ \HlJ \X2RJ \Hr 

and x~ is defined as 

fx; 

X, = — 
\Xi 

i.e. the charge-conjugate of (as required). 
Using the definition of the physical fields and A^, 



A^J \ sinUy^ cos9w J 



with 6^ being the weak mixing angle, we obtain from Eq. (jB.Sh for the left-handed parts 
of the x+i 
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coupling 


Feynman rule 


Xi r, X-i T 
Xi,RXj,R^tJ. 

xt xj 


-i e 6ij 7^ 



Table B.2: Feynman rules for X^X^ gauge boson couplings for an incoming Xj 



with 



^ JJ+-1 I 9 cos 6*^ 





' Cll 


Cl2 \ 




I Ci2 


C22 / 



The Cij are given by [75] : 



Cll 



C22 



gw 

cw 

gw 
cw 



2 3 1 

- 4 - 4 COS 2$^ 

2 3 1 

- 4 + 4 COS 2$. 



,Cl2 



gw 
Acw 



sin 2$_ 



For the coupling to the right-handed x"*", we obtain 



with 



-xf ( f f ] I'Z.xl, 

\ ^12 ^22 



I 

5 Cl2 





gw 
Cw 


„2 

Sw - 


3 

4 ~ 


1 

- cos 
4 


^22 = 


gw 
cw 


Sw 


3 

4 + 


1 

- cos 
4 



gw 
Acw 



sin 2$_ 



Similarly, the coupling of both left- and righthanded parts of are given by 

{xt\mt) = e 

with the term in £ mss m given by 

-^^je■X^l^A^l)t■ 
If we choose the x~ to be the particle, we have to substitute 

xl/Rl^X^xliR -Xr/l1^X^Xr/l- 
The Feynman rules for XX gauge-boson couplings are given in Table IB^ 
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B.4 Lepton-slepton and quark-squark chargino couplings 

We will now derive the general Feynman rules for the coupling of a sample SU{2) doublet 
to charginos. Feynman rules for MSSM particles can then be obtained by putting in the 
respective quantum numbers and the MSSM form of the superpotential. 
For a sample doublet 

the following couplings to charged gauginos and Higgsinos appear: 

• gauginos 

ig{ulW+dL + dlWuL) + h.c. 

(all SU(2) doublets) 

• Higgsinos 
terms of the form 

URVudhHi + URyudhHi +dB.yDULH^ + dnyoULH^ + h.c. 
Rewriting these terms using 

^ii^L^2 = mx2 , ^i-Pfl^2 = mxi, 

with ^' defined by [30] 
we end up with 

-g (i^uPRWdL + i^dPRW^'ui) + yu,ij ( i^d,jPRH^u*R^i + dL,j4^u,iPLH 



+yD,ij {i)u,jPRHd\i + ULjtpd^iPLH"^ + h.c. 



We use 



(Sl>'-'^(i)'(S><- 



— c ; ) 



with 



This yields 



\{-tW+)J \H+J \x. 
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(a) coupling to left-handed charginos 



(b) coupling to right-handed charginos 

{-gdhi^uiUi^Y + d*RyDi^u{Ur^Y)XiR 

(c) coupling to left-handed charge-conjugated charginos 

(d) coupling to right-handed charge-conjugated charginos 
(+ h.c. for all terms). 

They describe the following in- and outcoming charginos 

• incoming -|-, outgoing — 

(a), (6), (c)t, (d)t 

• incoming — , outgoing -|- 

(a)t, (6)t, (c), id) 

independent of the particle/ antiparticle choice. 
From the explicit for of the MSSM potential (lA.lip 

fussM = uVuQHu - dydQHd - eyeLHd + fJ-HuHd, 
we assign the following values to the sample matrices yD,yu used above: 





yu 


yo 


quarks 


yu 


yd 


leptons 





ye 


Higgs 









In the calculation, we assumed yu and yz) to be real (neglecting CP violating phases). In 
the high-energy limit, they are usually taken to be nonzero only for the third generations: 

yu oc J33, yn oc 533. 

The Feynman rules for an incoming xf are given in Table [ESI 
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(s) quarks 


(s)leptons 


Dol. 


Feynman rule 


(tbr) 




L 


iViiUt)* (a) 




(Ui L T ) 


R 


-ia(U^)* (h) 


OR) 


{^T Tr) 




^yb/T\^i2) \^) 


(bir) 


(t U-r r) 


R 


I Vk /^(U^)* (c) 


(duL), (scl), {btL) 


{Li Vi^L) 


L 


-^9iu+r (d) 


(btR) 




L 


^yt{ut2r (d) 



Table B.3: quark-squark and lepton-slepton couplings to xf • Vt/b/r are the (3,3) indices 
of the respective CKM matrices, k = e, /i, r. 

B.5 Chargino-Neutralino-Gauge boson and 
Chargino-Neutralino-Higgs couplings 

The part of Cmssm describing the chargino-neutralino-gauge boson couplings is given 
by 



fj, fJ, ^ 

from the gaugino sector and 
from the Higgsino sector. 

Rewriting this in four-component notation using the notations of the previous sections, 
we obtain H 

g (w^rW-W + -^{HdrW-PRH - Hu^^W^PlH)^ + h.c. 

The chargino-neutralino-Higgs coupling terms arise from the gauge part of the La- 
grangian 

-V2ig{'iijT''X''(t> - (f)*T''X''i;). 

Taking the quantum numbers of the Higgs doublets into account, they are given explicitly 
by 

^We use 
and 
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Hu terms: 



+ 



9 



V2 



+ 



V2 



+ h.c. 



Hd terms: 



MS' 
9 



+ {H-fW-H^, 



+ 



9 
\/2 



+ h.c. 



We can split this up into chargino-chargino-Higgs, neutralino-neutralino-Higgs, and 
chargino-neutralino-Higgs couplings. For actual calculations, we need to transform the 
Higgs bosons into Higgs mass eigenstates. We postpone this and here only sketch the 
coupling structure coming from the chargino/ neutralino sector. More details can be 
found in |142j . 

Chargino-Chargino-Higgs couplings 

The respective terms in the Lagrangian are given by 

ig{iH'j*W-H+ + (H'.yW+H-) + h.c. 
In four-component notation, this reads 

-g{WR{H^rHL + HR{H^^rWL) + h.c. 
Transforming this into the chargino eigenstates, we obtain 



9[Xlr 



Ur^{Hym^r + Ur^{H^^*m,r xlA+h.c. 



for the couplings of left- and right-handed charginos to Higgs bosons. Apart from the ro- 
tation into the Higgs mass eigenstates, these couplings agree with the ones given in [142] . 

Chargino-neutralino-Higgs couplings 

Chargino-neutralino-Higgs couplings are given by 

t({H+r{gW+Hl + j^W'H+ + ^BH+) 

+{H-ngW-H', - jf^'H- - ^BH,)) + h.c. 
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coupling 



Feynman rule 



xIrXIlH^ 



Table B.4: Feynman rules for an incoming x"* 



In four-component notation, we obtain 



^/2 



V2 



+ {Hd){-9Hd,LWR + ^WlHr + -^BlHr)) + h.c. 

Rotating this into the mass eigenstate basis, we obtain 

• Coupling of right-handed neutralinos to left-handed charginos 
^*^o / g' 



Coupling of left-handed neutralinos to right-handed charginos 



h.c. 



h.c. 



Apart from the rotation into Higgs mass eigenstates, this again agrees with the couplings 
as given in |142j . 

The Feynman rules for an incoming are given in Table [R!4l 
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C Point SPSla" 



700 

m [GeV] 
600 

500 _ 



SPSla' mass spectrum 



9 

9L : 



t2 



400 



X.4 



xt 



tl 



300 



200 



100 



II 
Ir 



T2 -^0 
X2 



^1 xS 



xt 



The mSugra point SPSla' used in this analysis is defined in [361 EZ]- It corresponds to 
a "typical" mSugra scenario. The respective mSugra parameters are 

Ml/2 = 250 GeV, Mq = 70 GeV, = -300 GeV, sign(^) = +1, tan/3 = 10. 

All sparticles masses in the point SPSla' are smaller than ITeV and therefore within 
the reach of future high energy colliders (LHC/ ILC). Branching ratios for decays can 
be found in [37] . 
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D Helicity eigenstates for massive fermions 



In the Dirac formalism, the helicity operator is given by 

for the projection of the spin along the Ck direction (see e.g. [143j ). Choosing the 
direction of flight along the z-axis, we can therefore define a spin up/ down projector by 

„ 1 ± acz 

^u/d 2 

As the helicity operator and the Hamiltonian commute, the Dirac spinors in the Dirac 
theory can be chosen to be eigenstates of the helicity operator and the Hamilton operator 
simultaneously; however, the helicity operator is not Lorentz-invariant0 and therefore 
frame-dependent. Samples of eigenstates can be found in [143] or [144] . 
The projector can be brought into a covariant form 

E(.) = 

for a spin vector s obeying 

P^s^ = 

(see [144j or |145j for more details). 
It can easily be shown that 

S(±s)S(±s) = S(±s) while S(±s)S(=Fs) = 0. 

The helicity projection operator as well as the energy projection operator then determine 
the eigenstates u±{p), Vzf{p) which solve the Dirac equation. For massless states (p^ = 0), 
the helicity projector reduces to the chirality projector Ph/l = • 
There are different approaches in the literature to deal with helicity states for massive 
fermions (see e.g. |146j for an original work or [147j and [127j for reviews on massive 
fermion treatment/ spinor techniques in general). However, we stick to the formalism 
introduced in [77]. Here, the authors construct u±{p), v±{p) which satisfy the Dirac 
equation and are eigenstates of the helicity operator. They then reformulate the 4- 
component S-matrix elements into spinor- matrix elements using Fierz identities, which 
they evaluate explicitly in dependence of the respective momenta in a given lab frame. 

^This only holds for massive particles; for massless particles, the helicity projector coincides with the 
chirality projector which is of course Lorentz-invariant. 
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Here, we just repeat the formulas used in the calculation of the helicity amplitudes and 
refer to the original work for further details. 

A four component matrix element can be reduced to a two component one using 
Here, 




Fierz identities give 

(W)a'7±(V'2)/3(V'3)7^/^±(V'4)5 = 0^1) a {i^ 4) 5 (4^1) j 0^2) (3 , 
(V'l)aO-±(V'2)/3(V'3)7^MT(^4)<5 = 2(V'l)a(V'2)/3 (V'3)7(V'4)<5 

where a, /3, 7, 5 = +, — . The product of two two-component spinors is given by 

{i^i)a{i^j)f3 = CiCjWaX,Wi3X^S{pi,Pj)x,X^, 

where 

Ck = l for {■lpk)r = {Uk)r, 
Ck = T for (V'fc)r = {Vk)T, 

w±{p) = Je± \^\. 



If S only depends on the momenta of the external particles, it reduces to the scalar 
quantity T: 

S{PiVj)\\j =T{piPj)x,x, = xUPi)X[^iPj) 
where x± ai^e helicity eigenstate spinors. Xij denote the helicities of the resulting massive 
particles in the given frame while a, /3 denote the chiralitie^. Tip iPj) are given as 

_ i\Pl \ +Pl,z){\P2 \ +P2,z) + {Pl,x -Wl,y){V2,x - W2,y) 



T{piP2)++ 
T{piP2)+~ 

T{piP2)~+ 
T{piP2)— 
If for example pi z - 



2\/|Pl|(IHI ^ P\,z)\P2\{\P2\ +P2,z) 
-{\pi \ +Pl,z)iP2,x -W2,y) + {\P2\ +P2,z){Pl,x 



iPl,y) 



2\/|Pl|(|Pl| + Pl,z)\P2\{\P2\ +P2,. 



-T*{piP2)- 
T*{piP2)+^ 



\Pl\ 



T{piP2)- 
T{piP2)^ 



P2,x + iP2,y 



^Wa\i should be read as WaxXi, i.e. w+ for a 



\/2|p5|(|p^| +P2,z)' 

Irnl +P2,z 
\/2|p^|(|H| +P2,z) 

X and W- for a X. W- = for massless particles. 
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Furthermore, it holds that 

which gives a complete determination of all possible T^^. 

The definition of the eigenstates given here always guarantees that the helicity eigenstates 
projected in an arbitrary frame correspond to eigenstates of the spin polarization in the 
rest frame of the respective particle; see [77] and [145j for further detail^. 



^You just have to identify the x± in IZZ] (3.19) with the A± in [145] (Eqn (2-50)). 
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E Generic diagrams contributing to NLO 
Chargino production 



In the following, we give the generic Feynman diagrams contributing to the process 
e'^e" — Xi Xi NLO. S denotes scalar, V vector particles, and F fermions. All 
diagrams including tadpoles are omitted. The output was generated using FeynArts. 
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e e ^ Xi %i 




e e ^ Xi %i 




e e ^ Xi %i 
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e e ^ xi Xi 




e e ^ xi Xi 
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F Soft and collinear approximation; ISR 
structure function 



F.l Soft Photon Factor 



The soft photon factor is given by Eq. ()3.6p . the integral appearing therein has been 
calculated in [87] and is given by 



d^k 2piPj apiPj 



|k|<AE^ 2a;fc {Pik){pjk) 



{api^ -p] 



1 , {aprf . 4Aii;2 

2 log ^ log ^ 



+ 



i 2 ^0 - U 

7 log -pT 

4 uo + |u| 



Li2 1 



Uo + |u| 



+ Li2[l 



Uq - |U| 



V 



u=api 



where 



and a defined by 



{ap, 



2{apio -PjoY 



a Pj — 2apiPj = , — > 0. 



Li2 are dilogarithms. The integral is regulated by the photon mass A. 



F.2 Hard-collinear approximation 
F.2.1 Finite mass effects in collinear radiation 

If photons are radiated off a charged particle under a very small angle 6, in the inte- 
gration over this region of phase space logarithms oc logm of the particle appear. The 
effects can be estimated by the collinear approximation, where in the squared matrix el- 
ement only terms oc -\ are kept and higher-order terms in are neglected. Work along 
these lines has been done in |148j : we just sketch the derivation and point to this work 
for further reference. We will deal with the case of unpolarized initial state particles and 
initial particles in a well-defined helicity eigenstate. 



106 



Unpolarized initial particle 



We start with the general description of an amphtude where one of the incoming particles 
radiates off a photon. Without the radiation, the amplitude is given by 

M = A{qi)u{p), 

where qi denote the momenta of the other particles involved in the process. If we include 
radiation, this cross section is modified to 

e 



M' 



2pk 



A{qi){f-l: + m)^{k)u{p), 



where k denotes the momentum of the photon and e its polarization vector. Squaring 
the amplitude and summing over the photon polarizations, we obtain 

For the photon polarization sum, we take advantage of the gauge freedom and choose 



qk 



where g is a general lightlike fourvector. We consider the g^^ and the ^ 

separately. 

Prom the g'^'^ term, we obtain 
\M\^ = 



part 



4 {pkf 

^2 



{pkf 



A{qi) (m-^ + {f- ^)m^ -kp{^ + m))A{qi), 



and for the second part 



A{pk) 



qk 



^kp — j(:pq + 2(m + p)pq 
qk 



-{m + i,)\A{qi). 



In the collinear approximation, we assume that 

\k \ . 



k 



P 



i.e. the angle 9 between A; and p is small. We then obtain 

nO 



k^ = 



X 



( P' \ 



V \i^\{i + o{e^)) J 



xp''il + 0{9)), 
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where 

^ = ^- 

Therefore, 

pk = pfk^ = xm2(l + 0(6*)). 

Now, we consider the terms pq and kq. In |148j . it is assumed that the angles a' between 
k and ~q, and f3' between 'p' and ~q are very large such that 

'k'^ = I'fc I I'^Ka + 0(0^)), 

'P^ = \-p\\^\{p + o{p^)), 

with 

a = - - a' , /? = 2 ~ 

and therefore 

^ = X (l + 0{a) + 0i(3) + 0{^ 



pq V (P° 

Note that this analysis also holds if ~q is nearly parallel to k and 'p' as for 



1 



%/2 V-fc 



Then, one has 



1 1 

- —p= \ k\ , p q = ^ p cos ^ 



and therefore 

/eg 2A;0 



We can then continue as before and replace 0{a) + 0{(3) by 0{9) in the error consid- 
erations. 

Now we look at the single contributions, where we only keep terms of C(m~^) and 
equally ignore higher-order terms in 0, a, [5. We obtain from the —g^^ term: 

\o{m-^){i + o{e))- 



{pkY (pk) 



2 



xp + o{e)o{ 



m 



-2\ 



kp kp 

UL _^ -o{m-^){i + o{e)) 

pk X 
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and therefore 

1 / / ,/x 9 , /,/ XN f (1 — x) x\ 

(ffl-^ -i- (lA — lt\m'^ — knlli- 4- mA\ = I 1 



i2 



For the ^'"'^ gfc'^'^^ part, one has 

1 

qk 



kqpk pk^ X 

pk kq 

and therefore 



]5A; \ gfc ^ J xpk ^ \m m? m? m 

Summing up everything, we arrive at 

pk ^^^^ (~ X pk 

as given in |148j . 



1^1' = ^x-^ MC'?^) (F-i) 



Initial particle in definite helicity state 

If the radiating particle is in a definite hehcity state, M' is given by 

A^' = A{q^) + ^(1 ± I^'Mp)^ 

where || /c and ^(1 ± 7^^) projects out the positive/ negative hehcity eigenstate. 
We parameterize s according to 

= pjP- - m^) + P^s^ (F.2) 
" m pq 

For the hehcity eigenstates, P± = 0, and we choose Py > and define positive/ negative 
hehcity by the it sign. 
We now have to calculate 
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with the same gauge choice as for the unpolarized case. We first consider the g^'' 
contribution and the and part of s^^, separately. Furthermore, we only calculate 
the part depending on (the part proportional to 1 in the helicity projector corresponds 
to the unpolarized case multiplied by ^). In the following, we use Z defined by 

\M'\^ = e^A{qi)ZA{qi). 

We consider the contributions coming from the combination of two terms in the polar- 
ization sum and two terms in s separately. In the following, the first term in the brackets 
denotes the term coming from the polarization sum ^ e^'^, the second the term stem- 
ming from s^. 



For the J- part, we obtain 



Z = -jp^ys (-2^-7^ [m^-kp) + m ■ f ■ -f^ - f ■ ^ ■ -f^ + A^i^ kp)) . 

We first look at the 0{w?) terms in the numerator; they are given by 

1 m? — kp ^ „ 1 c m? — kp , 1. 9^ .^.^n 

The leftover terms are of order ^0{m'~^). 
Prom the —"mj^ part, we obtain 

- 2^-7^ (m^ - m {kp)) + 2 m ^ • 7^ (kq) - 2 ■ 7^ {kq) - 2 m ^ • 7^ {pq) 
+ 2mf--f' [pq] + A-f' {kp) {pq)) . 

Here, the terms relevant for the approximation are given by 

7'(l--)(^-^) --^1' {x^-2x+l)f+{^ + -)0{m-^)0{e, a, 



2 pq 2 {pk)"^ ^ X 



\ kq ^ m J 
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rjij^g kf^q +^^^ part of the polarization sum is given by 

Z = /^i^J^Yqk i'^f'i'^^ ^^P^ + 2 m ^ • 1^ • 7^ (kp) -m^-p--f^ (hp) 

+ 2 ^ • 7^ (m^ - (kp)) (kp) - 2 ^ • 7^ (kq) - 2mp-^-'y^ (kq) 

+ 2m7^ (kp) (kq) +2^-7^ (kp) {kq) + 2m^ ^ • 7^ (pq) -m^-p--/^ (pq) 

+ mp-i(:-j^ (pq) +2^-j^ (kp) (pq) - • 7^ (kp) (pq)) . 

The 0{m'^) terms in the numerator are given by 



2{pk) 



pq 



1 / ^ , 2 f_ 1\\ (\ 1 



k'^q"+q>^k" mqP 
kq ' pq 

This part of Ztot is given by 
Z = -jj-^^{mp-^-t(:--f^-m^-i(:-p--f^+2^-p-j^{kp)-2m^--f^{kq) 
+ 2^-p-^^ {kq) - 2 ^ ^ • 7^ (pq)) . 
From this, we 

2 2-1 

■75(x-l)(#-|^) --^^5{x-l)p+-0{m-^)0{e). 



2{pk)^ ' ' ''-^ " 2(pA:)2 



Summing up all terms, we end up with 



2p/c \ j; pA; 

for the part proportional to Py as given in |148| . 

F.2.2 Helicity dependent structure functions 

In [HU [92l [93] , Eqs. (IF.lh and ()F.3p have been used to derive helicity-dependent structure 
functions which describe photon radiation integrated over the collinear region. Again, 



^Here we used that 
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we only sketch the result and refer to these works for further reference. 
Without photon radiation, the matrix element is given by 



and the square by 
With s'^ given by Eq. (|F.2p . we obtain 



and 

\M\^ = Aiq,) Q(l+/P||)^+0(„x))A(g, 
We now have to equate this with 

o2 



\M'\' = —A{q,){Ai + ^^A2j')A{q, 



where the Ai were derived in the previous sections. A short calculation then gives 



with 



B = )(2-2x + x'), 

X pk 

C = ^x\ (F.4) 
pk 



The next step is the approximation of ^ by 



pk 

1 2 



pk X (m2 + 4(p0)2 sin^D 

This approximation uses 

- p'k' (cos o-l^)+o{e')o (^) +o^^^^^ 



m 



Next, we have to approximate the integrals 
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in the collinear limit. Here, we use 

1 



m2 + 4(p0)2 sin^f m2 + (p0)2 ^2 

sine « + 0(03), 



^^,' m^ + {py{Aef \ ^ / (pO)2(A0)2 \ ^ m2 



^m2 m2 + (pO)2(A0)2y ^ m2 V V(^'°)MA0)2, 
As a next step, we approximate the integral over 

j dk^^\kl\ = J {pyx{l + 0{e))dx. 

Similarly, we assume 

j dTidT2dT3 6'''^\pi+p2-ki-k2-k3)F{pi,p2,ki,k2,h) = 

dk^ f / N 

^-J^^dn, J dT2dr3 6^*\p[+P2-k2-k3)Fip[,P2,k2,k3), 

when 

p'l = {Pi-h) = {I - X + 0(9)) pi. 
The differential phase space is given by 

d^h 



(27r)3 2A;0' 



The phase space integration approximation is exact for = 0. 

Taking all this into account, we then obtain for the integration over the factors B and 
C (Eq. dEl): 



|2 



/ \M^\'dx, 



I X 



2 2 

dVi ^. x^|.M~|^ = ^ / x\M-\'^dx. 
2 [pky 2 vr /„ . 

Here, denote the same/flipped radiating particle helicity amplitudes for the 2 
particle process and Xmin and Xmax are given by 

2AE^ 

Xmin = 1=-^ AE^: soft photon cut 

Vs 

(rn + m^Y 

Xmax = 1 kinematic limit for 3 body decay. 



min 

^max 
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If we want to bring this now in the form of a structure function 



a 



hard.coll i 



(s) = / (/+(x')<T+(x's) + /-(x')cT"(2;'s)) dx', 



we have to substitute 



X — s- x' = 1 — X 



x' 



and obtain 



^ ^ 27r X ^ ^ 



$(sx) / / vr? 

X 1 + C — 



The factor 

<I>(s) V~ ' V 

takes the difference of fluxes in the calculation of the total cross sections into account. 

F.2.3 Connection to leading log expressions 

In structure functions, collinear logarithms of the form 



In 



appear, where Q is the scale of the respective process. In these processes, Q is often 
taken as an upper limit for . We show that a change of variables from cos Q to pj_ 
exactly reproduces the results derived in section IF. 2. 21 and we can substitute 

in the structure functions. With k being the photon four-momentum 
kx_ = k^ sin0, 

COS^ = V^l^^=l-^(|52+0((^ 

d COS Q k I 



Furthermore, one has 



1 1 1 ( k\ \ 



pk pO/fcO 1 fci , 1 m2 V V(^°)V \W 
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Integrating 
we obtain 

i.e. exactly the result as given in the previous section. We therefore see that an upper 
limit 

exactly gives the usual structure function result integrated up to a scale 

F.3 ISR structure function 

F.3.1 Exact first order and soft solution 

In the following, we will consider the exact first order solution as well as the infinitely 
summed up soft solution for the electron structure function evolution equation (j3.13p 

1 dz 



with 

Pee{x) = , r? = log — 

\ 1 — X J TT 

(cf. Eqs. dam), ([311])). The +-distribution {h{x))+ is defined by 



(h(x)), = lim 



l-e 

e{l-x-e)h{x) -6{l-x) I h{y)dy 





(F.6) 



leading to 

" dxfix) (hix))^ = [ dx{f{x)- f{l))h{x). 

JO 

D^^ denotes the non-singlet part of the structure function [138| describing photon ra- 
diation only. In the derivation of the soft regime solution, we will closely follow [95j . 
The first order solution is given by 



£ - 



^iV5,a ^ lim (6{l-x) {- + 2 ln{e) ] + 9{1 - x - e 



o4 V V2 'V ' ' l-x 



l'\^^'''{x)dx = 1 + I (^2 ln(l - xo) + xo + ^xg^ (F. 
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The soft integrated 0{i]) contribution (|F.17p is negative because, when integrating from 
xq to 1, we do take all virtual, but not all real photon contributions into account. For 
xq = 0, the 0{ri) contribution vanishes and as expected, D^^''^{x) dx = 1. 

For the derivation of the solution in the soft regime, we perform a Mellin transform 
(cf . App ini) leading to 

= 1+1 r dpD''S<^^\p)P^^\ (F.8) 
4 Jo 

where 

= V'o(l) + ^^0(3) - MO - + 2). 

and tpoix) being the Digamma function defined in Appendix IG. 51 The formal solution 
to Eq. (jEH) is 

Z)^^(«)(r?) = expQr?P«)). 
Performing the inverse Mellin transform then leads to 

D'''iv,x) = ^exp (f(|-7)) £^JdCx-«exp(|(-^o(0 -^o(C + 2)) 

(F.9) 

Note that the exponential in front of the integral corresponds to the virtual part solution 
(modulo some constant terms which disappear when the ^ dependent digamma functions 
are subtracted). 
We then use 

MO + MC + 2) = 2 (^HO - E ^ + , (F.io) 

where i?2n denote the Bernoulli numbers (cf. Appendix lG.5|l . Only the logarithmic term 
of Eq. (jF.lOp is kept. This results in 



!.-(,,.). -exp^^(--,)jy^^^ .e^^exp^^(--,)^ ^^^^ 

(F.ll) 

Finally, we use 

- ln(y) = (1 - y) + 1(1 - y)2 + O ((1 _ yf) (F.12) 

and obtain 

=exp(|(2-.)) I <1^. (F.13) 
This is the solution to Eq. ()3.13p in the soft regime; it contains all orders in rj. 



116 



Finally, we can ask which splitting function would exactly lead to Eq. (jF.13p . We 
therefore look for a function /(^) such that (cf. Eq. ()F.9|1 ) 



J c—ioo ^ 



This is given by 

where for simplicity we used the Laplace instead of the Mellin transform. Therefore, we 
see that the approximation for y ~ 1 leads then to 

^-0(0 + ^-0(^ + 2) 2M0- 

This corresponds to 



1 - y 1 - y 

for the real photon part in the splitting function. In making this substitution, we omitted 
the term 

in the integration of the splitting function in the soft regime. We see that the (dominant) 
logarithmic behavior is exactly reproduced when integrating the modified structure func- 
tion in the soft region; the errors are 0(1 — y). 



F.3.2 Exponentiation 

The basic idea of exponentiation is to combine the 0{rf'),n oo emission of soft pho- 
tons, described in a leading logarithmic approach, with explicit finite order contribution 
in the hard-collinear regime into structure functions [1381 E3 EHl [961 1149[ I150[ . Then, 

• D{x) in the hard regime and 

• J^^ D{x) dx in the soft regime 

correspond to the exact solutions up to a certain order of r]. Comparisons of different 
approaches can be found in [75]. Following these lines, Skrzypek and Jadach obtained 
the following ISR structure function 
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/isr(2;) 



exp {-^VJE + I??) 



r(i + i) 



(l-x) 



+ 



16 

+ 4r 



-2(1 +2;)ln(l -x) 



21nx 



(1 + x) 



+-(l+x)lnx 



1 +x / 9 



vr 



X 5 
2 ~ 2 
1 



^ - ^ + - ln(l - x) + - ln2(l - x) - ^ In X ln(l - x) 



H In X 

16 



rLi2(l-x) 



+ 



- — Inx H In^ X Inxlnfl — x' 

8 12 2 ^ ' 



Here, 



1 + x^ 
2(1 -x) 

-^(1-x) nn(l-x) + ^J +^(5-3x)lnx + 0(7?^) (F.14) 



r, = ^ ln(% ) -1 . (F.15) 

This includes Eq. (IF. IIP as well as the exact solution for photon radiation up to 0{rf). 
We see that the definitions of rj in (|3.14|) and (|F.15p differ by a non-logarithmic factor; 
this corresponds to taking finite mass effects in the collinear approximation into account 
(cf. Appendix El)- 

Checking the requirements for exponentiation, we find that up to 0{r]), 



Wx) = |(l-x)i-i + -|(1 + 



x). 



(F.16) 



For the hard photon regime, we can expand this in r] and obtain 



fiSRix) 



1 



1 



1 -X 



(1 + x; 



1] 1 + X 



4 1-x 



This reproduces (IF.Sp . The term missing in the structure function expansion corre- 
sponds to the helicity flip contribution /~ which is subdominant. Integrating Eq. (|F.16p 
over the soft region yields 
"1 



fiSR dx 



(1 - xo 



XQ 



1 + 1^ 



V pi 1 2 

4 2""^°"^° 



1 + i 2 ln(l - xo) + xo + - xg + 0(r/2). 



(F.17) 

We see that this exactly corresponds to the result ()F.17p obtained from integrating the 
the 0{r]) splitting function Pee- We define the 0{a) part of Eq. ()F.17p fsoftjSR- 



/soft,iSR(^-E'7, ml) 



J f21n(l-xo)+xo + ^x2 ) . 



(F.18) 
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G Transformations and (Di)Gamma 
functions 



G.l The Fourier Transform 

The Fourier transform of a function f{x) is given by 

/oo 
-oo 

Its inverse is given by 

fix) =^ e^'"'fiP)dp, 
where the distribution of factors of ^== is purely conventional. Note that we used 



1 /""^ 
27r 7_oo 



G.2 The Laplace Transform 

The Laplace transform can be derived from the Fourier transform by taking the Fourier 
transform of 

g{x) = /(x)e-^^ 

where we assume that 

f{x) = for X < 
and additionahy making the change of variables 

s = c + ip. 
We end up with the Laplace transform 

/•oo 

C{s) = / e-"^/(x)dx 
Jo 

and its inverse 

1 rc+i oo 

fix) = —- / e'-Cis)ds. 

ZITT Jc—ioo 

Note that here the delta-function is given by 



1 



C+lOO 



S(x) = / e^^ds. 

2 ITT Jc—toc 

An exhaustive list of Laplace transforms and their inverses can be found in [151lll52"lll53j . 
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G.3 The Mellin Transform 



The Mellin transform can be derived from the Laplace transform by substituting 

y = e"^ ; g{y) = /(-ln(y)). 



We then obtain 



and 



UL 

s-1 



M{s) = / y'-'g{y)dy 







C+l oo 



9{y) = 7^ — / y 'M{s)ds. 



C—l CO 



Note that the upper limit of the Mellin transform integral UL depends on the value of 
f{x) for X < 0: 

1 for f{x) = for X < 
oo for f{x) 7^ for x < 0. 

In the first case, we can easily get the Mellin transform from the Laplace transform by 
performing the change of variables x = — ln(y). The delta function is here given by 



UL 



1 



c+i oo 



5(x — x') = / x ^ x'^ ^ds. 



c—i oo 



An overview on Mellin transformations can be found in |154j. 

G.4 The Euler Gamma function r(x) 

The Gamma function T{x) is defined by 

T{x) = / e~^e~^dt 



for X > 0. More useful representations, expansions, and relations to other functions can 
be found in [l55l [T56l[T57] . 

G.5 The Digamma function iIjq{x) 

The digamma function is formally defined by 

^o{x) = ln(r(x)) . 

A lot of useful representations, expansions, and relations can be equally found in the 
literature given in the last section. We here just quote a few properties useful in the ISR 
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derivation: 



Mx) = -7+ / \ 5/ (^>i)' (G.i) 







n— 1 ^ 

Tpoix + n) = '>po{x) + y~] — - 



fc = 



where 7 is the Euler constant. Eq. (|G.lh leads to 



t^-p^y^-Uy = -V.o(0-V'o(e + 2) + 2C, 
Jo 1-2/ 



where 



Another useful relation is given by 

B2n 



2£ ^2n£2n' 

n=l 



where B2n ai^e the Bernoulli numbers. 
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H References for Computer codes 

We list the main literature references for all computer codes mentioned in this work (in 
order of appearance) 

• WHIZARD [52] 

• FeynArts / FormCalc [HOl EH ESI IMl ES] 

• LoopTools |81j 

• CompHEP [nS] 

• MadGraph J20] 

• O'Mega [m] 

• CIRCE [122J 

• pythia ^5] 

• Sherpa [S3] 

• Isajet [54J 

• Form [1251 [T26] 
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